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ABSTRACT 


A  review  of  quantum  beats  in  beam-foil  spectro¬ 
scopy  is  presented.  Coherence  and  alignment  conditions 
are  defined  and  discussed  with  respect  to  the  obser¬ 
vation  of  zero-field,  Zeeman  and  Stark  beats.  The 
experimental  equipment  required  to  perform  quantum- 
beat  observation  is  described.  Methods  of  data  analysis, 
particularly  the  finite,  discrete  Fourier  transform  and 
direct  least-squares  approach,  are  developed  from  their 
fundamentals  and  applied  to  both  synthesized  and  actual 
data. 

Experiments  were  performed  to  observe  zero-field 
beats  in  boron  III  at  2066  A  and  Stark  beats  in  boron 
V  at  749  A.  Electric  field  calibration  was  carried 
out  by  the  observation  of  Stark  beats  in  hydrogen 
Lyman-a.  No  quantum  beats  were  conclusively  observed 
in  either  boron  experiment. 
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FOREWORD 


Beam-foil  Spectroscopy  at  the  University  of 
Alberta  was  begun  in  1967  by  Dr.  E.H.  Pinnington. 
Improvements  to  both  the  equipment  and  data  analysis 
occurred  steadily  through  1973  from  the  efforts  of 
staff  members,  a  former  Post  Doctoral  Fellow,  Dr. 

D.J.G.  Irwin,  and  former  students,  Drs.  C.C.  Lin  and 
A. E .  Livingston.  These  improvements  included  purchase 
of  a  new  vacuum  ultraviolet  spectrometer  and  photon¬ 
counting  detectors,  design  and  manufacture  of  several 
target  chambers,  laser  alignment  of  the  equipment, 
automatic  remote  data  acquisition  and  development  of 
an  extensive  data  analysis  program  for  lifetime  measure¬ 
ments  . 

The  author  joined  the  group  in  1973  and  parti¬ 
cipated  in  his  first  experiment  in  Spring  1974.  The 
purpose  of  the  experiment  was  the  first  attempt  at  the 
University  of  Alberta  to  observe  quantum  beats.  We  were 
fortunate  to  be  joined  by  Dr.  H.G.  Berry  of  the 
University  of  Chicago  and  Argonne  National  Laboratory 
(U.S.)  who  suggested,  during  the  course  of  the  experi¬ 
ment,  the  specific  attempt  to  observe  Stark  beats  in 
boron  V.  The  attempt  was  unsuccessful  and  Dr.  Berry 
invited  Dr.  Pinnington  and  myself  to  Argonne  National 
Laboratory  in  December  1974  to  try  again  with  higher 
beam  energy.  Again  we  failed  to  observe  quantum  beats 
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in  boron  but  did  observe  them  in  hydrogen  under 
special  experimental  conditions  which  had  only 
recently  been  reported  in  the  literature. 
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CHAPTER  I 


INTRODUCTION 


Beam-foil  spectroscopy  (BFS)  is  the  study  of 
radiation  emitted  by  a  beam  of  fast  moving  ions  after 
undergoing  excitation  by  passing  through  a  thin  target 
foil.  Pioneered  by  Stanley  Bashkin  (Ba  64)  at  the 
University  of  Arizona,  interest  soon  spread  to  other 
institutions,  particularly  following  the  first  inter¬ 
national  conference  on  BFS  in  1967.  Since  then,  two 
other  international  conferences  have  been  held  and  a 
fourth  is  scheduled  for  September  1975.  The  proceed¬ 
ings'*'  from  these  conferences  offer  an  excellent  review 
of  a  new  experimental  technique  in  contemporary  physics. 

The  beam-foil  light  source  is  characterized  by 
greater  advantages  and  fewer  drawbacks  than  conven¬ 
tional  spectroscopic  sources.  The  thin  foil,  typically 
3-15  yg/cm^  (^500  A)  of  carbon,  gives  a  sharply  defined 
spatial  origin  such  that  distance  "downstream"  from  the 
foil  corresponds  to  time  since  excitation.  The  low 
particle  density  eliminates  imprisonment  of  resonance 
radiation,  recombination  of  ions  and  electrons,  and 
stimulated  emission.  The  foil  excitation  easily  pro¬ 
duces  multiply-ionized  species  in  virtually  all  excited 

1  Ba  68,  MBB  70,  Ba  73. 
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states.  Magnetic  mass  resolution  effectively  isolates 

the  ion  of  interest  thereby  eliminating  impurities 

2  8  + 

except  in  cases  such  as  Si  which  can  be  contaminated 

with  N*  and  C0+  (typical  high  vacuum  residual  molecules). 

The  plethora  of  spectral  lines  which  are  usually 

observed  does  present  a  problem  with  line  blending. 

The  excitation  of  even  the  highest  energy  levels 

causes  cascading  to  the  level  of  interest  which  has 

proven  to  be  the  most  troublesome  drawback. 

One  obvious  application  of  such  a  light  source 

is  the  measurement  of  atomic  mean  lives.  Light  from 

a  short  segment  of  the  beam  (^0. 3-2.0  mm)  is  first 

viewed  as  close  to  the  foil  as  geometry  allows.  Pho- 

2 

tons  are  counted  for  a  preset  length  of  time  and  then 
the  foil  is  moved  a  calibrated  short  distance  and  the 
measurement  repeated.  By  successively  moving  down  the 
beam  the  measured  intensity  has  a  characteristic  expon¬ 
ential  decay  whose  decay  constant  yields  the  mean  life. 
Were  it  not  for  cascades  and  blends,  mean  lives  could 
be  as  reliably  measured  as  calibration  of  the  beam 
velocity  permitted.  Nevertheless  sophisticated  computer 
routines  and  diligent  analysis  by  the  experimenter  are 

2 

Actually  integrated  beam  current  in  a  Faraday  cup, 
or  light  normalization  is  preferable  to  time.  See 
Chapter  III. 
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able  to  provide  lifetime  values  for  almost  any  element 
in  several  stages  of  ionization.  In  many  instances 
these  measurements  can  be  converted  into  oscillator 
strengths  which  are  useful  guides  for  theoretical 
atomic  structure  calculations,  especially  for  heavier 
elements.  Another  recent  application  is  to  astro¬ 
physics  in  which  oscillator  strengths  are  needed  to 
derive  element  abundances  in  stars.  The  need  even 
extends  to  rare  earth  elements  whose  spectral  lines 
appear  to  a  remarkable  extent  in  photographic  plates 
of  chemically  anomalous  stars.  Such  stars  are  not  rare 
but  make  up  nearly  25%  of  the  main  sequence  including 
some  of  the  brightest  stars  in  the  sky. 

The  observation  of  light  modulations  along  the 
beam,  or  quantum  beats  (Chapter  II) ,  established  a 
different  family  of  measurements  that  could  be  made 
with  the  beam-foil  source.  These  include  fine-  and 
hyperf ine-structure  measurements,  Landd  g-values  and 
Lamb  shifts.  In  addition  to  such  frequency-related 
parameters,  quantum  beat  amplitudes  and  phases  may  be 
related  to  the  degree  of  non-isotropic  excitation  of 
the  atom  by  the  foil,  allowing  measurements  of  initial 
alignment  and  orientation. 

The  initial  plan  for  the  present  experiment  was 
to  observe  the  spectrum  of  boron  in  the  visible  and 
near  ultraviolet,  choose  a  selection  of  spectral  lines 
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which  showed  a  single,  gentle  exponential  decay  and 
look  for  zero-field  quantum  beats  between  hyperfine 
levels  since  that  has  yet  to  be  observed  in  boron. 
However,  the  failure  of  peripheral  equipment  forced 
us  to  convert  to  vacuum  ultraviolet  spectroscopy  and 
mean-life  measurements  while  the  equipment  was  being 
repaired.  We  looked  for  zero-field  beats  in  the 

O 

2066  A  radiation  of  boron  III  as  this  line  had  looked 
promising  in  a  previous  experiment  by  Dr.  Berry.  No 
beats  were  observed.  A  faint  line  of  boron  V  at  749  & 
had  also  been  detected  and,  being  a  hydrogenic  emission, 
a  Stark-effect  quantum  beat  measurement  could  have 
yielded  the  Lamb  shift  as  discussed  in  section  2.6. 

This  had  not  previously  been  measured  and  as  time 

3 

limitation  prevented  us  from  following  our  original 
plan,  we  concentrated  on  measuring  Stark  beats  in  the 
B  V  line. 

In  the  months  that  followed  we  attributed  failure 
to  observe  beats  to  a  poor  signal  to  noise  ratio.  This 
experiment  was  tried  again  at  Argonne  in  December  1974 
where  higher  beam  energy  allowed  us  to  collect  statis¬ 
tically  good  data.  In  fact,  although  the  statistics 
in  the  second  experiment  were  adequate,  our  analysis 


Accelerator  time  is  allocated  by  the  Radiation 


Research  Committee. 
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techniques  were  sufficiently  refined  for  us  to  know 
within  one  week  of  our  return  that  the  boron  V  data 
still  showed  no  beats.  Reasons  for  the  failure  of 
this  Lamb  shift  measurement  will  be  discussed  in 
Chapter  V.  As  explained  in  section  2.6  the  electric 
field  calibration  was  performed  by  examining  quantum 
beats  in  hydrogen  Ly-a.  This  offered  the  opportunity 
to  observe  "Eck  beats",  which  had  been  first  observed 
less  than  one  year  previously. 

Chapter  II  discusses  all  of  the  quantum  beat 
topics  mentioned  in  this  introduction.  It  is  pre¬ 
sented  as  a  general  review  and  does  not  employ  the 
mathematical  formalism  required  for  a  thorough  treat¬ 
ment  of  the  phenomena.  For  example,  the  use  of  the 
irreducible  tensor  operators  and  the  n-j  symbols  for 
angular  momenta  coupling  have  been  avoided.  The 
chapter  does,  however,  provide  a  framework  within 
which  the  varied  aspects  of  quantum  beats  may  be 
understood  and  a  good  starting  point  for  future,  more 
detailed  theoretical  investigations  of  this  complex 
subject. 

Chapter  III  is  a  brief  description  of  the 
experimental  equipment  used  for  this  project. 

Chapter  IV  describes  the  methods  developed  to 
analyze  the  data  and  represents  perhaps  80%  of  all 
the  time  spent  on  this  project.  As  is  indicated  in 
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the  text,  these  methods  need  to  be  extended  before 
they  are  generally  applicable.  However,  they  were 
quite  adequate  for  the  analysis  of  data  taken  for 
this  project.  Chapter  V  is  a  presentation  of  the 
results  of  these  experiments  and  a  discussion  of  the 
failures  encountered  along  the  way.  Criteria  for  the 
successful  observation  of  quantum  beats  in  future 
experiments  will  be  established. 
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CHAPTER  II 


REVIEW  OF  QUANTUM  BEATS 

Quantum  beats  are  modulations  of  the  radiation 
emitted  from  an  atom  after  it  has  been  excited  to  a 
coherent  superposition  of  upper  levels.  The  term 
"quantum  beats"  seems  to  have  been  applied  to  this 
phenomenon  when  first  observed  by  Aleksandrov  (A1  63) 
in  radiation  emitted  by  cadmium.  Following  the  theory 
of  Breit  (Br  33) ,  Aleksandrov  realized  that  the  band- 

■j  o 

width  of  the  3261  A  radiation  from  Cd  5s  XS  -  5p  p  was 

o  1 

broad  enough  to  encompass  the  Zeeman  splitting  of  the 

3 

upper  P^  term.  Consequently  he  irradiated  a  gas 

sample  of  cadmium  placed  in  a  magnetic  field  with  light 

from  a  cadmium  vapor  lamp  and  observed  quantum  beats  in 

the  radiation  from  the  sample. 

The  beam-foil  light  source  is  ideal  for  the 

observation  of  coherently  excited  states.  From  the 

uncertainty  principle  AE  equals  fi/At  where  At  for  beam- 

foil  is  the  time  the  atom  spends  in  the  foil,  typically 
-14 

10  sec.  This  gives  a  coherence  width  of  about 
.06  eV  (%500  cm  ^  or  ^15,000  GHz),  which  is  1,000  times 
greater  than  the  spatial  resolution  at  typical  beam 
velocities.  Coherence,  then,  is  guaranteed  by  the 
nature  of  the  excitation  for  all  quantum  beat  experi¬ 
ments  performed  so  far.  Further  requirements  for  the 
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observation  of  quantum  beats  will  be  discussed  in 
section  2.3. 

The  first  beam-foil  quantum  beat  experiment  was 
performed  by  Bashkin,  Bickel  and  Fink  (BBFW  65)  who 
oriented  the  entrance  slit  of  their  spectrometer  to 
view  lengthwise  along  the  beam.  They  observed  "bright 
spots"  on  the  photographic  plate  when  looking  at 
hydrogen  Lyman  radiation  for  n  =  3  to  7.  These  beats 
were  observed  in  a  very  weak  magnetic  field  of  8  gauss 
and  the  authors  attributed  the  beats  to  a  motional 
electric  field  producing  Stark  splitting  of  the  excited 
levels . 

Beam-foil  experiments  rapidly  became  more  sophis¬ 
ticated  with  the  introduction  of  photon  counting 
techniques  and  computer  analysis  of  data.  Since  vir¬ 
tually  all  experiments  centered  on  obtaining  lifetime 
data  it  is  interesting  to  note  that  the  first  theore¬ 
tical  treatment  of  quantum  beats  appeared  in  Macek ' s 
paper,  "Theory  of  Atomic  Lifetime  Measurements"  (Ma  70) 
in  which  the  author  pursued  the  necessary  conditions 
for  quantum  beats  in  order  to  determine  how  not  to 
observe  them.  Quantum  beats  were  considered  a  source 
of  error  for  lifetime  measurements!  But  Macek  recog¬ 
nized  the  interesting  consequences  of  his  work  which 
predicted  the  possibility  of  zero-field  beats 
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(section  2.4)  and  so  published  a  letter  to  that  effect 
(Ma  69)  prior  to  his  paper  on  lifetime  measurements. 
Andra’s  experiment  (An  70)  confirmed  Macek ' s  prediction 
and  represented  the  start  of  detailed  experiments  into 
the  quantum  beat  phenomenon. 

This  chapter  will  attempt  to  give  a  basic  under¬ 
standing  of  such  experiments  and  motivate  their 
inauguration  at  the  University  of  Alberta.  The  first 
section  describes  the  coherence  condition  and  derives 
quantum  beats  between  a  simple  two-level  system.  The 
third  section  treats  the  necessary  conditions  for 
quantum  beat  observation  using  the  density  matrix 
formalism  which  is  briefly  reviewed  in  section  2.2. 

The  last  three  sections  treat  zero-field  beats,  Zeeman 
beats  and  Stark  beats  respectively,  the  latter  repre¬ 
senting  the  type  of  experiment  performed  for  this  project. 

2 . 1  Coherence  Condition 

sudden  change 


- 1 - 

t  <  0  t=0  t  >  0 


Fig.  1.  Evolution  of  the  Hamiltonian 
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Basically,  the  quantum  mechanical  description  of 
quantum  beats  invokes  the  sudden  perturbation  approxi¬ 
mation  in  which  the  excited  atom  or  ion  is  considered 
to  be  in  stationary  states  of  a  Hamiltonian  from 
which  it  suddenly  enters  a  physical  region  described 
by  the  Hamiltonian  .  In  the  description  of  this 
latter  Hamiltonian  the  system  is  no  longer  in  a  sta¬ 
tionary  state  but  rather  a  coherent  superposition  of 
states  whose  time  evolution  gives  rise  to  the  beating 
phenomenon.  All  quantum  beat  experiments  following 
foil  excitation  can  be  classified  into  three  types 
according  to  the  Hamiltonians  H^  and  H2  used  to  des¬ 
cribe  the  system  (An  74) .  These  are: 

(1)  H^  describes  free  atoms  and  includes  external 
applied  fields. 

(2)  H^  describes  atoms  as  they  leave  the  foil  and  H2 
includes  external  applied  fields. 

(3)  H^  describes  atoms  as  they  leave  the  foil  and  H2 
describes  free  atoms. 

Type  (1)  is  applicable  to  systems  in  which  the 
external  field  region  is  some  distance  from  the  foil 
so  that  the  foil  excited  atom  has  decayed  to  a 
stationary,  excited  (perhaps  metastable)  state  from 
which  it  radiates  after  entering  the  field  region. 
Types  (2)  and  (3)  are  most  commonly  used  to  describe 
quantum  beat  experiments  in  which  the  field  region 
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(type  (2))  begins  at  the  foil.  These  are  called  Zeeman 
beats  for  magnetic  fields.  Stark  beats  for  electric 
fields,  and  appropriately  zero-field  beats  for  no 
external  field. 

The  simplest  quantum  mechanical  description  of 
a  system  which  predicts  an  oscillatory  modulation  of 
light  intensity  superposed  on  an  exponential  decay  is 
afforded  by  type  (1) .  A  hydrogen  atom  is  prepared  in 
the  | 2s>  state  and  suddenly  enters  an  electric  field. 
Following  Andra  (An  74)  the  wave  equations  without 
electric  field  can  be  written: 


2s (t) * 

|  2s>exp [- (ys/2 

+  ioo  )  t] 

o 

2p(t)>  = 

| 2p>exp [- (y  /2 

+  iojp)  t] 

with  |2s>  and  |2p>  solutions  of  the  time  independent 
Schroedinger  equation  and  the  time  evolution  factor 
including  a  phenomonological  damping  term  appears  in 
the  exponential.  In  like  fashion  the  solutions  in  the 
field  region  are  denoted  with  primes: 

l^2s,(t)>  =  I 2s 1 >exP I” (ys i /2  +  iws,)t] 

l^2p,(t)>  =  I 2P' >exPt- (ypi/2  +  iwp,)t]  . 

The  time  independent  solutions  can  be  expanded  accord¬ 
ing  to: 


. 


' 
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I  2s  •  >  =  a  |2s>  +  a  I  2p  > 

1  ss 1  ps  1  ^ 

(2.1) 

|2P'>  =  aspl2s>  +  appl2P>  * 

Assuming  orthonormality  of  the  field-free  states,  the 
expansion  coefficients  are  immediately  seen  to  be 

a  =  <2s|2s'>  a  =  <2p|2s’> 

ss  1  ps  ^  1 

a  =  <2s|2p'>  a  =  <2p|2p'>  . 

sp  1  pp  1 

Since  originally  the  system  was  prepared  in  the  1 2 s > 
state,  at  time  t  =  0  the  wavef unction  must  be  given 
by 

|2s>  =  |  2s ' ><2s 1  | 2s>  +  |  2p ' ><2p '  | 2s >  . 

At  later  times  the  wavefunction  is 

|^(t)>  =  a*  | 2s ' >exp [  -  ( y  , /2  +  iw _,)t] 

S  S  S  o 

+  asp I 2P ' >exp [~ (Yp i /2  +  iwpl)t]  . 

If  (2.1)  is  inserted  for  the  kets  on  the  right  this 


becomes 
|  ip  (t)  >  = 


+ 


(a 

(a 


* 


a 

ss  ss 

* 

a 

sp  sp 


2s>  +  assapsl2p>)eXp[-(^s,/2+  icV)t] 
2s>  +  aspappl2p>)exp["  (Yp,//2  +  i(V  }  t]  • 


The  reason  for  expressing  the  wave  function  in  this 
manner  becomes  apparent  from  the  expression  for  the 
intensity  of  spontaneous  emission: 


I  (t)  =  Const  |  <ls  |  dipole  operator  |  i/j  (t)  >  | 

From  this  it  can  be  seen  that  transitions  can  occur 
only  through  the  | 2p>  admixture.  Finally  this  yields 

I (t)  =  A  exp (~Y  , t)  +  B  exp (-y  , t) 

s  p 

+  C  exp[-(ys,/2  +  yp  ,/2 )  t]  cos  (u>s  , -u)p  ,  )  t 

which  shows  a  modulation  whose  frequency  yields  the 
separation  of  the  |2s’>,  | 2p ' >  levels.  If  the  electric 

field  were  calibrated  and  measurements  were  made  at 
several  field  values,  an  extrapolation  to  zero  field 
would  yield  a  measurement  of  the  Lamb  shift  (section  2.6) 
The  basic  concept  of  quantum  beats  is  that  a 
sudden  change  of  basis  allows  the  observation  of  a 
sinusoidal  time  evolution  of  the  quantum  states.  The 
above  discussion  exemplifies  this  but  is  too  cumbersome 
for  the  treatment  of  more  complex  systems.  Before 
tackling  the  problem  with  greater  generality  a  brief 
review  of  the  density  matrix  will  be  given. 

2 . 2  Density  Matrix 

Perhaps  the  most  thorough  description  of  the 
density  matrix  was  given  by  Fano  (Fa  57)  in  the  Reviews 
of  Modern  Physics.  Some  of  the  fundamental  aspects 
will  be  given  here  following  his  review. 


■ 


' 
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Consider  a  state  vector  \ p  whose  eigenfunction 
expansion  is 

ib  =  T  c  u 
Y  L  n  n 
n 

with  the  un  orthonormal.  The  mean  value  of  any  operator 
Q  is  given  by 


< 


Q>  =  y  Q_,_  c*  c 

4  n '  n  n '  n 
n  1  n 


If  a  nonpure  state  is  represented  by  an  incoherent  super¬ 
position  of  pure  states,  ^ ^ , to  each  pure  state  corres¬ 
ponds  a  mean  value  <Q>^,  each  with  statistical  weight 
p  ^  ,  and  the  mean  value  of  Q  for  the  incoherent  super¬ 
position  is  given  by  the  grand  average 

<Q>  =  l  p(l)<Q>.  =  l  Q  ,  l  p(l)  c(;J-)*c^l) 

4  ^  1  n’n  4  ^  n’  n 

l  nn '  l 

If  one  relates  the  second  term  to  the  density  matrix 
then 

V  (i)  (i)*  (i) 

P  i  =/P  C  |C 
rnn'  4  ^  n'  n 

l 


and 


<Q>  J  ^n'npnn'  ^,^p^n'n 
n '  n  n ' 


or 


<Q>  =  Tr(Qp) 


(2.2) 


where  Tr  indicates  the  trace  of  the  matrix.  Equation 
(2.2)  is  taken  as  the  definition  of  p. 
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Several  useful  properties  follow  from  the  above 
discussion: 


(1) 

(2) 

(3) 

(4) 

(5) 


Pn'n  =  pnn,?  the  densitY  matrix  is  Hermitian 
Tr(p)  =  1 

2  0  in  any  representation 
Tr  ( p2 )  <  1 

2 

For  pure  states  (p  )  =  p 

mn  mn 


Time  evolution  is  given  by 


=  4-  [H,P] 

at  ^ 


-mt/fi  iHt/fi 

p  (t)  =  e  '  p  (0 )  e  ' 

<Q>t  =  Tr [Q p  ( t ) ]  =  Tr [Q  (t) p  (0 ) ]  . 


2 . 3  Alignment  Condition 

The  intensity  of  light  emitted  in  the  transition 
from  a  number  of  upper  states,  denoted  i ,  to  a  number 
of  lower  states,  f,  is  given  by  (FM  73)^ 

I(t) =  C  l  | <f |e* .D | i (t) >|2  (2.3) 

f,i  ~  ~ 

with  e  the  unit  real  or  complex  polarization  vector,  D 


This  is  a  common  but  not  simple  result  and  many  older 
texts  offer  no  derivation  whatever;  see  Mi  70,  Chap.  4 
and  Lo  73,  Chap.  8. 


■ 


. 
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the  operator  (usually  dipole)  inducing  transitions  and 
C  is  a  proportionality  factor  including  those  terms 
needed  to  make  I  an  energy  flux.  It  could,  for  example, 
include  an  integration  over  e  if  unpolarized  light  were 
being  observed.  Inherent  in  (2.3)  are  the  assumptions 
that  any  splittings  within  the  upper  or  lower  states  are 
small  compared  to  the  gross  structure  involved  in  the 
transitions,  and  that  all  included  transitions  fall  well 
within  the  bandwidth  of  the  detection  system  (AD  73) . 

If  (2.3)  is  expanded  in  an  eigenbasis  of  the  Hamiltonian 
H2  of  section  2.1  it  becomes 

I(t)=C  Y  <f  |e*.D  I  a><«  |i  (t)  ><i  (t)  I  3><3  I  e.D  I  f  >.  (2.4) 

f /I 

OLr  3 

If  £|i(t)><i(t) |  is  recognized  as  the  initial  density 
1 

matrix,  p(t),  and  ^  |  f  > < f |  =  p_  the  projector  onto  final 

f  f 

states  (2.4)  becomes 

I(t)=C  I  <3|  (e.D)p,(e*.D)  |a><a|p(t)  |  3>  .  (2.5) 

a, 3  ~  ~  f  ~  ~ 

From  the  previous  section  the  evolution  of  the  density 
matrix  is 

p(t)  =  U(t)  p(0)U+(t) 

with 

U(t)  =  exp  (-  £  Ht) 

H|a>  =  ft(wa  -  j  ya) |a>  . 


and 


. 


. 


' 


The  presence  of  a  phenomenological  damping  term,  y  , 

CL 

requires  mathematical  justification  which  shall  not  be 
given  here  but  see  BG  74.  y^  can  be  taken  as  the 
inverse  of  the  lifetime  of  the  state  |a>.  The  inten¬ 
sity  then  becomes 

I(t)  =  C  7  <  3|(e.D)p_(e*.D)  |a><a|p(0)  |  3> 

3^  ~  II  <%*  ~ 

. 

x  exp  [-i  (a)a-C0g)  t]  exp  [-^(y^+y^)  t]  .  (2.6) 

This  is  the  fundamental  expression  for  quantum 
beats.  It  shows  that  beats  are  possible  so  long  as 
(MG  74) 

^<3 | e.D | f ><f | e* .D | a>  ^  0  (2.7) 

and 

<a | p  (0)  | 3>  ?  0  .  (2.8) 

The  first  is  the  condition  that  each  of  the  upper 
states  must  be  able  to  decay  to  the  same  final  state 
via  the  same  transition  operator  D.  The  energy 

<v 

structure  of  final  states  does  not  affect  the  beat 
pattern.  The  second  is  the  condition  that  the  initial 
density  matrix  in  the  eigenbasis  of  E^  must  have  off 
diagonal  elements.  This  is  known  as  the  alignment 


condition. 


'  ^ 

•; 


■ 
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Certain  assumptions  about  the  initial  density 

matrix  can  significantly  simplify  the  problem. 

Because  the  foil  interaction  takes  place  in  such  a 

-14 

short  time  ( < 1 0  s)  the  density  matrix  may  be  written 
in  an  eigenbasis  of  composed  of  uncoupled  atomic 
states  |LMLSMgIMI>  if  the  atom  obeys  LS  coupling. 
Furthermore  in  this  representation  p(0)  is  assumed  to 
be  spin  independent  so  that 

<LMlSMsIM].  |  p  (0)  |  L  ’M^SMglM-J. 

<lml | p (0) |L'M^> 

(2S+1)  (21+1)  * 

If  the  foil  is  perpendicular  to  the  beam,  axial  symmetry 
demands  that  p(0)  be  invariant  for  any  rotation  about 
z-axis  chosen  as  the  beam  axis.  It  must  also  be  invar¬ 
iant  under  reflection  in  any  plane  containing  the  beam. 
These  symmetry  conditions  lead  to  (Ma  70) 

<lml|p(0) |l'm^>=  6m  m,<l|ml|  |p(0) |l' |m^|> 

L  L 

from  which  it  can  be  seen  that  p(0)  in  the  represen¬ 
tation  is  diagonal  for  a  given  L.  The  alignment  con¬ 
dition  demands  that  these  diagonal  elements  must  not 
all  be  the  same.  This  can  be  shown  by  writing 

<a|p(0)|B>  =  1  <a|i>a.<i|8> 

i 

If  were  the  same  for  all 


with  o ^  =  <i|p(0)|i>. 


' 


■ 
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diagonal  elements 


<a | p (0 )  | 3>  =  a  J<a  |i  ><i | 3>  =  a<a|3>  =  0 

i 

since  the  sets  | a>  and  |i>  are  taken  to  span  the  same 
subspace.  But  <a|p(0)  | $>  ^  0  by  (2.8)  so  not  all  the 
ct  are  the  same. 

As  a  specific  example  of  this  alignment  condi- 

2 

tion  consider  the  n  =  2  P  state  of  hydrogen  in  zero 
field  ignoring  hyperfine  structure.  Then 

<  Ik  MlMs  |  p  (0  )  \  lk  M^Mg  >  = 


ml 

1 

0 

-1 

k 

-3* 

k 

-k 

k 

°1/2 

0 

0 

0 

0 

0 

1 

-k 

0 

o±/2 

0 

0 

0 

0 

k 

0 

0 

V2 

0 

0 

0 

0 

-k 

0 

0 

0 

a  /2 
o 

0 

0 

-1 

k 

0 

0 

0 

0 

V2 

0 

-k 

0 

0 

0 

0 

0 

V2 

H2 

for  this 

case  is 

the 

free  atom 

Hamiltonian 

with  LS 

coupling.  The  transformation  to  the  J,  MT  represen- 

u 

tation  is  accomplished  with  Clebsch-Gordon  coefficients: 


. 
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<1  JMT  |  p(0)  |1  3^  j'm!>  =  V  <  1 3g  JM  T  1 1 3g  Mt  M_  > 

MM  J  L  S 

L  S 

X<1  ^MlMs  |  p  (0  )  1 13s;  MlMs><1  3gMLMs  |l  3g  J'Mj>  . 

All  of  these  coefficients  may  be  generated  from  the 
simple  table  (CS  63)  . 

< j x 3^  |  j1k  jm> 


After  evaluating  the  coefficients  and  performing  the 
indicated  sum,  the  transformed  density  matrix  may  be 


written 


. 


' 
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<1  h  JMt  I  p (0 )  |l  h  j'm'  >  = 

J  J 


This  shows  by  example  that  off-diagonal  elements  occur 
only  if  gq  ^  o^.  Note  too  that  off-diagonal  elements 
only  occur  for  like  values  of  MT. 

To  summarize  this  section,  equations  (2.7)  and 
(2.8)  represent  necessary  conditions  for  quantum  beats. 
Each  of  the  upper  states  must  couple  to  the  same  ground 
state  by  the  same  transition  operator.  The  initial 
density  matrix  is  solely  a  function  of  the  beam-foil 
interaction.  After  considering  the  short  excitation 
time  in  the  foil,  spin  independence  is  postulated  and 
that  plus  axial  symmetry  and  reflection  symmetry  leads 
to  a  density  matrix  diagonal  in  the  uncoupled  atomic 
representation.  That  the  diagonal  elements  are  not 
all  equal  is  the  alignment  condition  verified  by  the 


. 


' 

■ 
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existence  of  zero-field  quantum  beats.  The  existence 
of  off-diagonal  matrix  elements  in  the  basis  states 
of  ^  also  follows  from  the  sudden  perturbation  con¬ 
dition  (coherence)  which  for  zero  field  allows  the 
transformation  to  be  made  via  Clebsch-Gordon  coeffi¬ 
cients.  The  density  matrix  elements,  a^,  are  referred 
to  as  excitation  cross  sections  since  they  are  depen¬ 
dent  only  on  the  foil  interaction. 

2 . 4  Zero-field  Beats 

The  existence  of  zero-field  quantum  beats  was 
first  suggested  by  Joseph  Macek  (Ma  69)  and  first 
observed  by  H.J.  AndrS  (An  70)  who  followed  the  sugges¬ 
tion.  It  should  be  apparent  from  the  previous  section 
that  zero-field  beats,  with  an  assumption  of  no  cas¬ 
cades,  are  amenable  to  a  thorough  theoretical  treatment 
since  transformations  required  to  evaluate  matrix 
elements  involve  only  well  known  coupling  coefficients. 
The  complications  which  arise  are  purely  spatial  and 
thorough  treatments  have  been  given  in  the  literature, 
most  notable  being  FM  73,  El  73  and  BG  74. 

Equation  (2.5)  could  be  written 

I (t)  «  tr [Mp (t) ] 


with  M  the  monitoring  operator  (CP  66)  proportional  to 


. 


. 


, 

.. 


4* 

(e.D)  (e.D)  .  If  M  and  p  were  both  expanded  into 

spherical  tensor  operators  the  various  moments  of  p 
could  be  related  to  combinations  of  different  pro¬ 
ducts  of  J.  For  example  alignment  with  axial  symmetry 

•  2  2 
is  the  same  as  a  nonzero  expectation  value  of  3J  -  J  . 

z 

If  the  expectation  value  of  J  were  nonzero  the  system 
would  have  initial  orientation  (possible,  for  example, 
if  axial  symmetry  were  broken) .  The  formalism  allows 
one  to  relate  various  combinations  of  the  polarization 
of  the  emitted  light  (Stokes  parameters)  to  the  various 
excitation  cross  sections  permitted  by  the  foil  inter¬ 
action. 

To  obtain  all  the  parameters  of  interest  experi¬ 
mentally,  the  emitted  radiation  should  be  viewed 
through  a  quarter-wave  plate  and  linear  polarizer. 

The  polarizer  is  set  parallel,  perpendicular  and  at 
45°  to  the  beam  and  this  latter  measurement  repeated 
after  rotating  the  quarter-wave  plate  45°.  This 
entire  set  of  measurements  is  then  repeated  with  the 
optical  axis  forming  a  different  angle  to  the  beam. 

Just  such  a  meticulous  experiment  has  recently  been 
performed  by  Berry  et  al.  (BCES  74)  in  conjunction  with 
a  tilted  foil  experiment.  The  results  showed  a  clearly 
measurable  effect  on  the  polarization  as  the  foil 
angle  was  changed  justifying  in  principle  the  need 
for  a  more  powerful  formalism  than  the  one  of  the 


- 


t' 


. 

. 


. 
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previous  section. 

It  should  be  emphasized  that  with  zero-field 
beats  the  total  intensity  is  not  modulated  and  the 
experiment  is  best  performed  with  a  polarizer.  This 
was  confirmed  by  Burns  and  Hancock  (BH  73)  who  viewed 
the  2s  S  -  3p  P  3889  A  radiation  from  He  I .  By 
giving  careful  consideration  to  the  foil  position, 
they  concluded  that  light  polarized  parallel  to  the 
beam  was  a  pure  cosine,  twice  as  intense,  and  180° 
out  of  phase  with  light  polarized  perpendicular  to  the 
beam.  With  the  polarizer  set  to  54.8°  relative  to  the 
beam  axis  such  that  the  parallel  and  perpendicular 
components  were  combined  in  the  ratio  1:2,  the  signal 
was  proportional  to  the  total  intensity  and  was  not 
modulated.  These  results  are  in  strict  agreement  with 
theory  which  postulates  a  spin  independent  foil  exci¬ 
tation  giving  strong  support  to  Macek's  analysis. 

The  modulations,  then,  for  zero-field  are  traced  to 
a  reversible  exchange  of  angular  momentum  between 
internal  degrees  of  freedom  (FM  73) . 

2 . 5  Zeeman  Beats 

The  Zeeman  effect  is  the  splitting  of  an  atomic 
level  into  its  component  magnetic  sublevels  by  an 
external  uniform  magnetic  field  B.  The  Hamiltonian 
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of  an  atom  in  such  a  field  is  given  by 


H  =  H  +  -r—  (L  +  g  S )  .  B  , 
a  h  ~  ^s  ~  ~ 

with  the  Bohr  magneton  and  gg  the  factor  by  which 

the  electron  magnetic  dipole  moment  differs  from  that 

of  a  classical  particle  with  angular  momentum  S.  If 

the  field  direction  is  taken  as  the  z-axis  the  diagonal 

matrix  elements  of  L  and  S  are  (Mi  70) 

z  z 

<LSJMj|Lz|LSJMj>  =  flMj  S(S+1) 

|  r>  Itot*,  .  _  S(S+1)  +  J(J+1)  -  L(L+1) 

<LSJMT  S  LSJMt>  =  tlM-,  - ^  \ - 

J  1  z 1  J  J  2J  (J+l) 


First  order  perturbation  then  gives  the  splittings  of 
the  energy  level 

<aLSJMT  |H|aLSJM_>  =  e  (aLSJ)  +  g_y_M_B  (2.9) 

J  J  J  O  J 


with 


i(i  +  gs)+  (gs-  i) 


S  (S+l)  -  L  (L+l) 
2 J (J+l) 


Since  g  s  2  within  1  part  in  one  thousand 
s 

„  .  3  ,  S  (S+l )  -  L  (L+l )  . 

gJ  ~  2  2J (J+l)  ' 


written  in  this  form  gT  is  called  the  Land£  g-factor. 

J 

The  energy  separation  between  adjacent  magnetic  sub- 


levels  from  (2.9)  is 


' 

• 

. 


. 
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ft(UL  =  gjyBB  ' 

with  co  =  2 tt  times  the  Larmor  frequency  which  is  the 
classical  precession  frequency  of  a  magnetic  dipole, 
with  moment  gTy_.,  in  an  external  field. 

The  necessary  conditions  to  observe  quantum  beats 
are  satisfied  if  the  beam  is  subjected  to  a  uniform 
transverse  magnetic  field  which  may  extend  right  to  the 
foil.  Due  to  the  longitudinal  alignment  along  the  beam 
direction  a  magnetic  field  perpendicular  to  the  beam 
induces  coherence  between  states  with  AM  =  2.  The 

u 

light  intensity  polarized  parallel  to  the  beam  (if  the 
magnetic  field  is  at  right  angles  to  both  the  beam  and 
the  axis  of  observation)  may  then  be  written  (MG  74) 

I  (t)  =  I  (0)  [  1  +  A  cos  [2ooT  t)  ]  exp  (-t/x)  .  (2.10) 

A  measure  of  the  beat  frequency  gives  a  measure  of  the 
Land£  g-f actor,  a  parameter  of  fundamental  importance 
for  determining  the  validity  of  the  coupling  scheme 
employed  for  the  description  of  atomic  states. 

The  usual  experimental  technique  employed  with 
Zeeman  beats  is  to  observe  the  beam  at  a  fixed  dis¬ 
tance  from  the  foil  and  to  linearly  change  the  mag¬ 
netic  field  thereby  eliminating  the  exponential  decay. 
In  order  to  maintain  a  small  magnetic  field  and  still 
observe  a  sufficient  number  of  oscillations,  the  chosen 


■ 


. 
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detection  position  should  be  as  far  downstream  as 
possible  though  not  so  far  that  the  signal  intensity 
is  lost  in  background  noise. 

Though  Zeeman  beats  offer  the  possibility  of 
measuring  Land£  g-factors  for  ionization  states 
unattainable  by  other  means,  there  is  one  severe 
drawback  to  this  method  being  applied  to  any  atom. 

To  produce  observable  splittings  small  magnetic  fields 
(^100  gauss)  are  employed  but  for  many  atoms  the 
splittings  so  induced  are  comparable  in  magnitude 
to  the  hyperfine  structure  caused  by  the  nuclear  spin. 
All  experiments  performed  so  far  which  use  equation 
(2.10)  have  measured  samples  with  1=0.  When  1^0 
the  theoretical  description  becomes  very  complicated 
(GCBL  73)  and  the  experimental  results  of  Gaillard 
et  al.  differed  greatly  from  their  theoretical  predic¬ 
tions  for  3  out  of  4  measured  levels. 

2 . 6  Stark  Beats 

The  application  of  a  static  electric  field 
produces  mixing  of  states  with  different  ^-values 
for  hydrogenic  atoms.  States  with  the  same  n  and  j 
quantum  numbers  are  energetically  degenerate  even  in 
the  Dirac  theory.  However  in  a  now  famous  experiment 
Lamb  (La  52)  measured  an  energy  difference  between 


1 
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2  2 

the  2s  S2  and  2p  Pa  levels  in  hydrogen  which  has 

'i  '2. 

come  to  be  called  the  Lamb  shift.  Quantum  electro¬ 
dynamics  (QED)  is  able  to  predict  an  energy  shift 
from  the  Dirac  value  for  all  fine  structure  levels 
based  on  an  interaction  between  the  electron  and  the 
quantized  radiation  field.  Measurements  of  the  Lamb 
shift,  usually  performed  with  resonance  techniques, 
are  among  the  most  precise  in  all  of  physics.  The 
amazing  agreement  between  these  measurements  and  the 
predictions  of  QED  is  perhaps  the  greatest  achievement 
of  quantum  mechanics. 

The  observation  of  Stark  beats  with  beam-foil 
spectroscopy  allows  direct  measurement  of  the  Lamb 
shift  as  a  zero  field  interpolation  of  energy  split¬ 
tings  obtained  at  several  electric  field  values.  Though 
the  accuracy  cannot  hope  to  rival  that  of  resonance 
techniques,  it  is  possible,  at  least  in  principle,  to 
measure  the  Lamb  shift  between  levels  which  are  not 
metastable  in  a  variety  of  hydrogenic  ions  unavailable 
before  beam-foil. 

The  major  aspect  of  our  experimental  attempt 

2 

was  to  measure  the  Lamb  shift  between  the  4p  P2/2~ 

2 

4d  D3/2  levels  in  boron  V  by  looking  at  the  n  =  4  to  3 
transition  at  749  A.  Fig.  2.  is  a  partial  energy 
level  diagram  for  boron  V  excluding  hyperfine  structure 
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but  including  both  radiative  corrections  and  corrections 
for  finite  nuclear  charge  distribution  and  finite 
nuclear  mass.  The  calculation  was  performed  using  the 
equations  compiled  by  Garcia  and  Mack  (GM  65) . 

When  a  hydrogen  atom  is  subjected  to  an  electric 
field  the  linear  Stark  effect  produces  splittings  in 
the  fine  structure  levels  corresponding  to  different 
values  of  m|j|»  The  seven  fine  structure  levels  of 
n  =  4  split  into  16  levels  which  yield  a  maximum  number 
of  120  different  frequencies!  However  many  of  these 
are  nearly  identical  and  most  do  not  fall  into  the 
experimental  resolution.  A  program  to  evaluate  the 
energy  levels  for  given  electric  field  had  previously 


been  written  here  following  the  procedure  of  Luders 
(Lii  51)  .  The  output  was  used  to  determine  which  fre¬ 
quencies  might  be  observable  and  to  plot  an  extrapola 
tion  curve  for  those  that  were.  The  results  will 
be  discussed  in  Chapter  V.  The  possibility  that  non¬ 
linear  effects  might  be  important  at  fields  as  high 
as  1000  volts/cm  was  not  considered. 

In  order  to  calibrate  the  external  electric 
field  which  was  axially  applied  (see  Chapter  III) 
the  experiment  was  performed  with  hydrogen  and  the 

O 

Lyman-a  transition  at  1216  A  was  observed.  Quantum 
beats  were  detected  at  two  different  electric  field 
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values  and  the  resulting  calibration  agreed  favorably 
with  geometrical  considerations. 

The  field  calibration  with  Lyman-a  offered  an 
experimental  opportunity  to  check  a  different  aspect 
of  Stark  beats.  Though  the  frequency  patterns  in 
quantum  beats  can  be  quite  complicated,  in  principle 
they  can  be  predicted  theoretically.  Measurements  of 
frequencies  simply  confirm  theoretical  models  or  in 
the  case  of  Zeeman  beats  shed  light  on  the  proper 
coupling  scheme  to  use.  Quantum  beat  amplitude  and 
phases  however  cannot  yet  be  theoretically  calculated 
since  they  depend  in  an  unknown  manner  on  the  beam- 
foil  interaction.  Most  current  QB  research  is  directed 
toward  understanding  this  interaction.  Recently  Eck 
(Ec  73)  proposed  an  experiment  for  hydrogen  in  which 
data  is  collected  with  an  electric  field  first  parallel 
then  anti-parallel  to  the  beam  axis.  The  experiment 
was  performed  by  Gaupp  et  al.  (GAM  74)  who  interpreted 
their  results  to  mean  the  electron  cloud,  in  zero  field, 
is  displaced  in  front  of  the  proton  as  the  atom  leaves 
the  foil.  It  was  a  simple  matter  for  us  to  perform 
this  experiment  while  at  Argonne  but  unfortunately  time 
only  permitted  one  set  of  data  to  be  taken  with  reversed 
electric  field. 

The  discussion  in  sections  2.3  and  2.4  centered 
on  quantum  beats  between  levels  of  a  single  atomic 
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term  such  as  P.  For  all  nonhydrogenic  atoms  the 
energy  separation  between  terms  is  too  large  to 
spatially  resolve  though  not  necessarily  too  large 
for  coherence  to  be  present.  But  for  hydrogen, 
different  terms  are  essentially  degenerate  and 
though  symmetry  arguments  predict  a  density  matrix 
initially  diagonal  in  and  Mg  for  given  L,  there 
is  no  reason,  as  Macek  pointed  out  (Ma  70)  why  off- 
diagonal  elements  for  different  L  but  like  and 
could  not  exist.  Two  types  of  coherence  may  be 
distinguished,  that  involving  states  with  equal 
parity  such  as  S-D  coherence,  and  that  involving 
states  with  odd  parity  such  as  S-P  coherence.  The 
former  was  actually  observed  in  hydrogen  Lyman- 3  under 
zero-field  conditions  by  Burns  and  Hancock  (BH  71) . 

The  latter  can  only  be  observed  using  the  Stark  effect 
which  mixes  terms  of  opposite  parity. 

The  relevant  portion  of  the  density  matrix  for 
n  =  2,  with  Mc  suppressed  since  M2  =  M  a  ,  is 

a 
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Eck  claims  that  the  beat  intensity  is  proportional 
to  (Ec  73) 

2 

I(t)  «  &  [i(a  +  2 q  )  -  a  ]  cos  tot 
to  3  p0  p,  s 


Vco 

o  V 

— o —  Re  a  cos  tot  +  -  Im  a  sin  ut  , 

id2  A  sp  U/3  sp 


where  V  is  the  Stark  matrix  element  (in  angular  frequency 

units)  coupling  the  Sx  m. =  \  and  P1  m.  =  ^  states,  and 
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co  =  (toQ  +  4V  )  2  with  (oQ  the  Lamb  shift.  V  is  simply  a 
constant  times  the  electric  field  intensity.  If  data 
is  collected  with  the  field  first  parallel  then  anti¬ 
parallel  to  the  beam  and  the  beats  are  added  and  sub¬ 
tracted  one  obtains: 
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If  the  diagonal  elements  of  p  are  considered  known 

from  previous  measurements  (DAWB  72,  AD  73),  the 

amplitude  and  phase  of  the  difference  signal  allows 

evaluation  of  a  .  The  calculation  with  results  will 

sp 

be  presented  in  Chapter  V. 
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CHAPTER  III 


EXPERIMENTAL  EQUIPMENT  AND  TECHNIQUES 

The  data  acquisition  phase  of  the  quantum  beat 
experiments  used  the  same  experimental  arrangement  as 
has  previously  been  used  to  collect  lifetime  data. 

The  original  beam-foil  set  up  is  described  in  Li  71. 

The  current  arrangement  contains  extensive  modifica¬ 
tions  from  the  early  work  to  be  found  in  Li  74  and 
IL  75.  The  description  given  here  will  be  a  brief 
overview  of  all  the  equipment  including  mention  of  the 
on-line  computer,  Texas  Instrument  980  A,  used  for  the 
first  time.  The  description  is  that  of  the  University 
of  Alberta  Radiation  Research  Laboratory  and  while 
facilities  at  Argonne  National  Laboratory  (US)  are 
much  more  extensive,  an  identical  spectrometer  and 
nearly  identical  target  chamber  were  used  there,  so 
no  separate  description  of  that  facility  will  be  given. 
Differences  where  applicable  will  however  be  mentioned. 

Fig.  3  is  a  diagram  of  the  experimental 
arrangement  used  here.  The  particle  accelerator  is  a 
2  MV  Van  de  Graaff  using  radio-frequency  excitation  of 
a  chosen  source  gas.  (At  Argonne  it  is  a  4  MV  Dynami- 
tron) .  Two  separate  source  gas  containers,  thermo- 
mechanically  controlled,  allow  two  different  samples 
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ION  BENDING 


Fig. 3  Experimental  Arrangement 


Fig. 4  Beam  Length  Acceptance  Geometry 
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to  be  studied  without  opening  the  machine.  Boron 
trifluoride  and  hydrogen  were  used  for  this  experi¬ 
ment.  The  r-f  excitation  normally  produces  singly 
ionized  species  and  H+,  H+ ,  H*,  B+,  BF+,  and  BF* 
were  all  observed.  Mass  separation  is  accomplished 
by  an  electromagnet  which  deflects  the  beam  11°.  With 
a  full  strength  field  of  5  kilogauss  a  singly-charged 
particle  of  200  a.m.u.  is  so  deflected.  The  two 
stable  isotopes  of  boron,  ^B  and  ^^B,  were  resolved 
and  the  more  abundant  ^'*"B  was  used  for  studies  at 
749  A. 

After  deflection  in  the  magnet,  the  mass-selected 

ion  beam  enters  the  target  chamber,  which  is  under 

reasonably  high  vacuum  (10  torr) .  There  it  passes 

2 

through  a  thin  carbon  foil  (5  yg/cm  )  and  terminates 
in  a  Faraday  cup.  Beam  currents  of  1  yA  are  typical 
but  for  boron  at  Alberta  the  maximum  current  was  0.25 
yA. 


The  Target  Chambers 

Three  different  target  chambers,  each  with  its 
own  advantages  and  disadvantages,  were  considered  for 
this  experiment;  two  at  Alberta  and  one  at  Argonne. 

Of  the  two  at  Alberta  one  is  a  focused  system,  described 
in  Li  71,  employing  a  fixed  parabolic  mirror  and  a  small 
rotatable  plane  mirror.  The  target  foil  remains 
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stationary  and  the  beam  is  scanned  by  rotating  the 
plane  mirror.  Six  foils  are  mounted  on  a  wheel 
allowing  targets  to  be  changed  without  breaking  the 
vacuum.  Presently  the  mirror  coatings  are  such  that 
this  chamber  cannot  be  used  in  the  far  vacuum  ultra¬ 
violet  region  of  the  spectrum  but  plans  to  have  the 
mirrors  platinum  coated  are  being  considered.  Though 
the  chamber  was  mounted  and  ready  for  use  at  the  time 
of  the  experiment,  the  failure  of  peripheral  equipment 
forced  a  change-over  to  the  lifetime  chamber  used  for 
mean  life  studies  in  the  vacuum  ultraviolet. 

The  second  chamber  at  Alberta,  used  for  all 
recent  mean  life  studies,  is  non-focusing  and  contains 
a  target  holder  mounted  on  a  rotating  screw.  The  beam 
axis  is  as  close  to  the  spectrometer  entrance  slit  as 
possible  (6  cm)  and  the  beam  is  scanned  by  physically 
moving  the  target  upstream.  The  target  moves  under 
control  of  a  small  stepping  motor  which  is  geared  to 
the  mounting  screw.  Control  of  the  motor  is  handled 
remotely  by  a  somewhat  sophisticated  electronic  cir¬ 
cuit,  which  allows  fully  automatic  or  manual  control. 
Physical  position  of  the  gears  is  visually  monitored 
with  a  digital  turn  count  indicator.  Variable  step 
sizes,  as  small  as  0.1  mm,  are  possible.  This  target 
chamber  also  allows  foils  to  be  changed  under  vacuum; 
as  many  as  15  foils  are  available.  Though  an  electric 
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field  can  be  applied  to  the  beam,  the  foil  changing 
mechanism  prevents  the  field  from  being  uniform. 
Experiments  involving  the  use  of  an  electric  field 
must  therefore  allow  for  field  calibration. 

The  target  chamber  at  Argonne,  designed  with 
this  experiment  in  mind,  has  good  electric  field 
geometry.  The  foil  is  mounted  across  the  end  of  a 
small  hollow  cyclinder  which  is  then  inserted  into 
a  large  aluminum  plate  (about  1  cm  thick) .  This  large 
plate  is  then  mounted  to  a  screw  drive  which  is  driven 
in  the  same  manner  as  the  one  at  Alberta.  A  second 
large  but  thin  plate  is  mounted  with  insulators  2  cm 
downstream  from  the  foil  plate  and  the  beam  passes 
through  a  small  hole  in  the  plate.  Opposite  potentials 
up  to  5,000  volts  were  applied  to  these  plates  to 
observe  Stark-effect  quantum  beats.  Unfortunately 
this  chamber  could  accommodate  only  one  foil  at  a  time. 

Spectrometer  and  Detectors 

The  spectrometer  employed  at  both  Alberta  and 
Argonne  is  a  McPherson  Model  225  one-meter  mono¬ 
chromator.  The  grating  used  is  a  Bausch  and  Lomb 
aluminum  replica  (600  rulings/mm)  overcoated  with 
magnesium  fluoride  and  with  platinum.  The  grating 
is  concave,  96  mm  wide  and  56  mm  high,  blazed  at  900  A 
in  first  order.  The  entrance  and  exit  slits  are  at  the 
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radius  of  curvature  of  the  grating.  Since  there  is 
no  further  focusing,  the  length  of  beam  from  which 
light  can  reach  the  detector  (call  it  £)  is  solely  a 
function  of  slit  width,  s,  grating  width,  w,  and 
distance  of  the  beam  axis  from  the  entrance  slit,  which 
as  mentioned  previously  is  6  cm.  Since  the  distance 
from  slit  to  grating  is  one  meter,  elementary  geometry 
yields  the  formula  for  £  (Fig. 4)  : 

t  =  s  +  10--+  ”>■  .  (3.1) 

100 

By  masking  the  grating,  w  can  be  set  to  any  desired 
value.  At  Alberta  two  masks  for  grating  widths  of 
10  and  20  mm  are  frequently  used  and  the  option  of 
changing  mask  size  without  opening  the  spectrometer  is 
available.  The  slit  width  is  externally  adjusted  to 
any  value  with  a  micrometer  screw.  Of  course,  improve¬ 
ment  in  resolution  by  masking  the  grating  and  narrowing 
the  slit  is  paid  for  by  a  reduction  in  light  intensity. 
For  this  experiment  two  combinations  were  used.  The 
best  resolution  was  with  s  =  0.6  mm  and  w  =  10  mm 
giving  £  =  1.2  mm.  Some  of  the  boron  data  were  taken 
with  s  =  0.8  mm,  w  =  20  mm  giving  £  =  2.0  mm. 

The  detector  used  below  1100  A  is  a  Bendix  4219 
Spiraltron  electron  multiplier  with  no  window.  The 
dark  count  for  this  detector  is  less  than  0.1/sec. 
Between  1050  A  and  2000  A  an  EMR  542  G  photomultiplier 
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with  lithium  fluoride  window  is  used.  Its  dark  count 
is  typically  0.4/sec. 

Data  Collection 

Before  any  data  collection  actually  begins, 
several  spectral  scans  at  different  beam  energies  in 
the  wavelength  region  of  interest  are  recorded.  The 
spectrometer  grating  can  be  turned  by  hand  or  with  a 
built  in  motor  with  several  available  speeds.  The 
detector  output  is  connected  to  the  pen  of  a  strip- 
chart  recorder  which  produces  the  spectral  scan. 

Several  well-known  lines  are  always  available  for  cali¬ 
bration  of  the  output.  Different  ionization  states  are 
identified  by  noting  whether  the  intensity  of  a  line 
increases  or  decreases  with  changing  beam  energy.  Since 
Boron  V  is  hydrogenic,  the  n  =  4  to  n  =  3  transition 
energy  can  be  easily  calculated  to  be  at  749  A.  This 
line  does  not  begin  to  make  an  appearance  until  a  beam 
energy  of  1.8  MeV  is  reached.  Once  a  line  has  been 
selected  for  study,  the  spectrometer  is  left  untouched 
and  data  collection  begins  from  a  remote  control  center. 

Four  parameters  specify  a  given  datum  during  any 
one  run  of  a  beam-foil  experiment;  they  are, 

(1)  foil  position,  which  represents  the  time  since 

ion  excitation  at  the  foil 


. 
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(2)  photon  counts,  the  light  intensity  integrated 
over  the  counting  time 

(3)  counting  time,  the  length  of  time  during  which 
the  detector  signal  is  accepted  by  the  counter 

(4)  integrated  beam  intensity,  the  beam  current, 
digitized  and  summed  over  the  counting  time. 

The  stepping  motor  circuit  can  be  triggered  by  either 
a  preset  time  (time  normalization)  or  a  preset  current 
(current  normalization) .  To  compensate  for  fluctuations 
in  the  beam,  current  normalization  is  usually  chosen. 
Three  parameters  are  then  recorded:  foil  position, 
photon  counts,  and  counting  time.  The  latter  value  is 
used  to  subtract  the  dark  count  from  the  data.  Record¬ 
ing  of  the  foil  position  is  most  useful  for  lifetime 
experiments  in  which  small  steps  are  used  initially, 
while  larger  ones  are  used  near  the  end  of  the  run. 

For  quantum  beats  equal  step  sizes  are  used  throughout 
the  run  so  this  parameter  becomes  redundant. 

Two  methods  have  been  used  at  Alberta  for  record¬ 
ing  the  data,  but  recent  developments  with  the  TI  980  A 
system  promise  to  improve  these  methods.  The  first 
employs  a  Teletype  ASR  33  which  simultaneously  punches 
a  paper  tape  and  produces  a  teletype  print-out.  After 
the  experiment  the  paper  tape  is  converted  to  computer 
cards  which  are  edited  to  remove  false  starts  and 
other  mistakes,  and  then  run  through  a  data  condensation 
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routine  which  produces  cards  suitable  for  analysis 
by  the  program  HOMER,  as  is  discussed  at  the  beginning 
of  the  next  chapter.  The  second  method  uses  a  Fabritek 
multichannel  analyzer,  capable  of  storing  256  numbers 
in  each  of  four  quarters.  Data  in  the  Fabritek,  at  the 
time  of  the  last  experiment,  could  be  read  by  the 
TI  980  A,  although  there  was  no  way  to  produce  an 
output  suitable  for  the  IBM  360/67.  Now,  however, 
this  shortcoming  is  soon  to  be  removed  and  our  data 
recording  methods  will  be  changed  accordingly. 

On-line  Computer 

At  Alberta  a  Texas  Instrument  980  A  computer  with 
a  Silent  700  terminal  was  used  during  data  acquisition 
to  perform  a  preliminary  analysis.  At  the  time  of  the 
experiment  only  a  magnetic  tape  drive  was  available  as 
peripheral  equipment  and  core  memory  was  limited  to 
16  thousand  bytes.  Many  routine  computer  events  (such 
as  compiling  a  FORTRAN  program)  were  performed  by 
reading  from  and  writing  to  two  cassette  tapes  on  the 
Silent  700.  A  100  statement  program  required  40 
minutes  to  compile.  Since  then  a  disk  drive  has  been 
added,  memory  has  been  increased  to  24  K  and  an  X-Y 
plotter  is  operational.  Compiling  time  is  now  less 
than  one  minute  for  the  previously  mentioned  program. 
The  Silent  700  can  also  now  be  used  as  a  terminal  to 


- 


' 


. 


, 


43 


the  IBM  360/67  and  soon  the  TI  980  A  will  be  able  to 
read  and  write  MTS  files.  This  final  improvement  on 
the  system  will  greatly  enhance  the  Radiation  Research 
facility . 

At  the  Argonne  National  Laboratory  an  ASI  2100 
computer  and  peripheral  equipment  has  been  functional 
for  several  years.  All  stages  of  data  acquisition, 
including  triggering  of  the  stepping  motor,  were 
handled  by  the  computer.  Data  was  collected  on  mag¬ 
netic  tape,  printed  with  a  line  printer  and  punched 
onto  IBM  cards  before  the  experiment  was  concluded. 

The  computer,  however,  was  totally  inadequate  for  any 
attempt  at  preliminary  data  analysis  and  this  contri¬ 
buted  to  the  failure  of  the  experiment.  Had  it  been 
obvious  that  there  were  no  quantum  beats  on  data  taken 
at  an  optimum  field  voltage  with  a  beam  acceptance 
length  of  2  mm,  the  grating  would  have  been  masked  and 
the  slit  width  changed  to  achieve  maximum  resolution. 
This  change  would  have  more  than  doubled  the  counting 
time  thereby  reducing  the  number  of  field  voltages  at 
which  data  could  be  taken.  In  fact,  this  change  was 
made,  after  data  had  been  taken  for  14  hours  at  five 
optimum  field  voltages.  It  then  took  11  hours  to 
collect  data  at  the  highest  and  lowest  field  values, 
the  former  having  too  high  a  beat  frequency  and  the 
latter  too  small  a  beat  amplitude  to  resolve  the  beats. 
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Preliminary  data  analysis  for  a  quantum  beat  experiment 
is  most  desirable. 


CHAPTER  IV 


METHODS  OF  DATA  ANALYSIS 


Two  different  techniques  are  discussed  in  this 
chapter  for  analyzing  beam-foil  quantum  beat  data. 

The  first  is  a  Fourier  transform  method  and  the  second 
is  a  direct  least-squares  fit.  The  analysis  was  per¬ 
formed  on  an  IBM  360/67  computer  and  the  programs, 
written  in  FORTRAN  IV,  are  listed  in  the  appendix. 

Since  the  quantum  beats  appear  as  small  modulations 
superposed  on  a  multi-exponential  background,  both  of 
the  above  techniques  enlist  the  aid  of  a  lengthy  program, 
called  HOMER.  HOMER  was  developed  here  at  the  University 
of  Alberta  between  1970  and  1972  to  extract  lifetimes 
from  beam-foil  data  and  is  briefly  described  in  section 
4.1.  The  data  analysis  proceeds  in  two  stages.  The 
first  is  the  evaluation  of  the  background  using  HOMER 
and  the  second  is  the  evaluation  of  the  frequency,  am¬ 
plitude  and  phase  of  any  quantum  beat  using  both  the 
Fourier  transform  and  least-squares  methods. 

Two  short  auxiliary  programs  have  also  been  used 
to  test  the  success  of  the  analysis  methods.  The  first 
is  one  to  synthesize  data.  This  uses  applicable 
experimental  parameters,  up  to  five  background  para¬ 
meters  (double  exponential  plus  constant) ,  and  up  to 
ten  modulation  parameters  (five  frequencies,  five 
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amplitudes  and  alternately  cosine  and  sine  phase) .  The 
second  is  a  program  to  impart  "random"  scatter  to  any 
data  sample,  regardless  of  whether  the  data  is  synthetic 
or  actual.  The  program  is  referred  to  as  the  routine, 
FUZZ,  and  has  proven  especially  useful  in  determining 
signal  to  noise  characteristics  of  a  data  sample.  Of 
course,  no  scatter  which  is  generated  within  a  program 
can  be  truly  random  and  this  routine  is  limited  to  using 
one  of  five  random  numbers  with  which  it  generates 
noise.  Which  of  the  five  numbers  is  used  is  selected 
by  the  user  prior  to  running  the  program.  The  magnitude 
of  the  noise  is  determined  by  the  square  root  of  the 
number  to  be  altered.  FUZZ  was  lifted  in  its  entirety 
from  HOMER  where  it  appears  as  an  optional  subroutine, 
and  is  listed  in  the  appendix  for  the  interested  reader. 

All  of  the  output  appearing  in  this  chapter  is 
artificial,  generated  by  the  above  mentioned  programs. 
Whether  or  not  the  data  was  "fuzzed"  will  be  made  clear 
in  the  text. 

4 . 1  Finite  Discrete  Fourier  Transform 

A  useful  comparison  between  the  usual  and  finite 
discrete  Fourier  transform  is  given  in  GS  66.  The 

usual  definition  for  the  continuous  case  is 

00 

/• 


F  (k) 


—  CO 


f  (x)  exp(2irikx)  dx 


,1  — 


•  • 
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while  for  the  discrete  case 

N 

F  (m)  =  l  f  (n)  exp  (2  7rinm/N)  .  (4.1) 

n=l 

N  is  the  total  number  of  data  points,  f (n)  the  datum 
at  position  n,  while  F (m)  is  the  transform.  Both  n 
and  m  are  integers  running  from  1  to  N.  To  obtain  the 
inverse  transforms  one  only  needs  the  orthogonality 

conditions  which  for  the  continuous  case  is 

00 

■ 

exp  [27rik  (x-x  '  )  ]  dk  =  6(x-x') 

* 

—  00 

while  for  the  discrete  case  it  is 
N 

\  exp  [2tt  im  (n-n '  )  /N]  =  N6M(n-n')  (4.2) 

m=l  N 

where  6..  is  the  Kronecker  delta  function  with  its 
N 

argument  considered  modulo  N;  that  is,  6^(pN)  =  1  for 
integer  p,  otherwise  6^  =  0.  The  inverse  transforms 
can  easily  be  obtained  in  the  usual  manner  yielding  for 
the  discrete  case 

T  N 

f(n)  =  [  F  (m)  exp  (-2tt  inm/N)  .  (4.3) 

m=l 

Note  that  for  the  continuous  case  the  usual  factor  of 
2tt  does  not  appear  outside  the  integral  because  of  its 
inclusion  in  the  exponential.  This  has  been  done 
deliberately  so  that  k  would  equal  1/A,  the  wavenumber, 
commonly  used  in  spectroscopy. 


. 


-V 


To  relate  the  integers  n  and  m  to  physical 
parameters  the  arguments  of  the  exponentials  must 


be  related  by 


(4.4) 


with  t  the  time  since  foil  excitation,  v  the  frequency 
k  the  wavenumber  and  x  =  ct;  note  that  x  is  not  the 
physical  position  downstream  from  the  foil.  This 
physical  position  (call  it  x')  is  related  to  the  beam 
velocity,  v£,  and  the  step  size,  d,  by  x'  = v^t  =  nd. 
Therefore 


(4.5) 


x 


n 


c 


where  S  is  called  the  scale  factor  and  is  constant 


c 


over  any  data  set.  Typically  0.5  cm  <  S  <50  cm. 

c 

With  this  value  for  x  eq.  (4.4)  gives  k  =  m/N S  .  The 

c 

frequency  spectrum  is  discreet  as  expected  starting  at 

k(low)  =  1/N  S  .  The  upper  limit  would  appear  to  be 

m  =  N  but  for  the  fact  that  the  second  half  of  the 

transform  is,  except  for  phase,  the  same  as  the  first 

half.  The  upper  limit  thus  occurs  at  m  =  N/2.  Then 

k(high)  =  1/2S  which  corresponds  to  each  data  point 

o 

being  alternately  above  and  below  the  background. 
Trustworthy  results  for  experiments  with  fewer  than 
100  data  points  are  not  obtained  close  to  either  limit 


; 


- 
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The  modulations  under  consideration  are  superposed 
on  a  background  which  must  be  removed  before  the  trans¬ 
form  analysis  can  be  performed.  Calling  this  background 
Bk(t),  the  data  can  be  decomposed  as 

Data  =  R ( t )  Bk  (t)  +  Bk  (t) 

with  R(t)  the  normalized  residuals  (beat).  The  only 
requirement  for  Bk  is  that  it  be  a  least-squares 
approximation  to  the  data  which  ignores  the  modulation. 
For  this  analysis  Bk  was  assumed  to  be  a  sum  of  expon¬ 
entials  and  its  determination  was  accomplished  by  the 
HOMER  routine  (IL  74)  which  can  fit  a  maximum  of  six 
parameters  in  the  following  form: 

Bk(t)  =  A1  exp  (-t/x^)  +  A2  exp  (-t/x2 )  +  A^  expC-t/x^) 

where  t,  as  before,  is  the  time  since  foil  excitation 
and  A^ ,  x^,  A 2,  A^ ,  x^  are  determined  by  the  pro¬ 

gram.  The  optimum  number  of  parameters  (between  2  and 
6)  is  that  which  minimizes  the  reduced  chi-squared. 
Criteria  for  goodness-of-f it  as  well  as  a  clear  over¬ 
all  discussion  of  the  multi-exponential  data  reduction 
problem  can  be  found  in  Li  74. 

For  computer  analysis,  sine  and  cosine  functions 
are  substituted  for  the  complex  exponentials  in  (4.1). 
Each  are  added  and  then  squared  to  produce  a  graph  of 


■ 


:  ' 


F(k)  F*(k).  The  individual  sine  squared  and  cosine 
squared  values  are  printed  but  not  graphed.  They  are 
used  to  determine  the  phase  of  the  beat  as  explained 
in  a  later  paragraph.  The  individual  transforms 
actually  used  for  data  analysis  are 

N 

F  (k)  =  V  R(n)  cos  (2-rrnS  k)  (4.6a) 

o  ^ «  c 

n=l 

N 

F  (k)  =  Y  R(n)  sin(27rnS  k)  (4.6b) 

S  n=l  C 

with  F  F*  =  F^(k)  =  F^  +  F^.  When  evaluating  this  at 

c  s 

a  selected  set  of  values  for  k  it  is  not  necessary  to 
restrict  the  plot  to  only  N/2  points.  Intermediate 
values  may  be  evaluated  to  produce  a  smooth  transform 
but  the  resolution  limit  is  set  by  the  orthogonality 
condition.  If  a  peak  is  observed  at  some  value  k 
corresponding  to  an  integer  m,  then  the  amplitude  will 
be  zero  at  m  ± 1.  The  half  width  of  the  peak  is  then 
Am  =  1  or  Ak  =  l/NSc. 

Fig.  5  shows  a  semilog  plot  of  synthetic  data 
with  a  cosine  beat  at  2.25  GHz  (5%  amplitude)  and  a 
background  function 

B (t ) =  7000  exp(-t/0.1  ns) +  3500  exp(-t/0.4  ns) . 

O 

This  background  resembles  that  of  B  V  at  749  A.  The 
scale  factor  corresponding  to  typical  experimental 


. 
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Counts 
x  10~  3 


Fig. 5  Simulation  of  Boron  V  Data 


. 


conditions  was  chosen  to  be  1.0  cm  and  40  points  were 
generated.  The  resolution  is  Ak  =  0.025  cm  ^  which 
is  also  the  starting  frequency  corresponding  to  m  =  l. 
Fig.  7  shows  the  Fourier  transform  in  the  low  fre¬ 
quency  part  of  the  spectrum.  The  maximum  amplitude 
occurs  at  m  =  3  while  the  amplitude  is  zero  at  m  =  2 
(1.5  GHz)  and  m  =  4  (3.0  GHz). 


4 . 2  Amplitude  and  Phase  Determination 

For  sine  and  cosine  functions  orthogonality  may 
be  written: 


y  cos(2irkx)  cos  (2-rrk ' x)  =  j  6N(k-k') 
x=l 


l  sin(2irkx)  sin  (27rk '  x)  =  ^  6N(k-k')  (4.7) 

x=l 


N 

£  sin(27rkx)  cos  (27rk*  1  x)  =  0  . 
x=l 

Properly  the  sum  is  over  n  not  x,  and  the  argument  of 
the  Kronecker  delta  is  m  not  k.  Assuming  the  form  of 
the  modulation  to  be  A  cos  (27rk '  x-0 )  ,  (4.6a)  gives 

N 

F  (k)  =  I  A  cos  (27Tk'x-9)  cos(27Tkx) 
c  u , 

x=l 

N 

=  y  A  cos  0  cos(2Trk'x)  cos  (2Trkx) 
x=l 

N 

+  y  A  sin0  sinC27Tk'x)  cos  (2iTkx)  . 
x=l 
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SUMSQD 

SUMSIN 

SUMCOS 

FREQU. 

ENERGY 

(GHZ) 

(CM-1 ) 
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0.0250 
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-0. 4804E-01 

0.900 
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-0 . 3550E- 01 
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0.0325 
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0 . 1 034E-02 

1  .050 
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1  .200 
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- 0 . 9664E- 0 1 

0. 1310E+00 

1.275 

0.0425 

0 . 1739E-01 

-0. 4132E-01 

0.1 252E+00 

1.350 

0.0450 

0 . 62  78E- 02 

0.3105E-02 

0 . 7917E-01 

1  .425 

0.0475 

0. 1523E-03 

0.  1228E-0  1 

0. 1194E-02 

1  .500 

0.0500 

0. 8556E-02 

-0.3202E-01 

-0. 8678E-01 

1.575 

0.0525 
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-0. 1671E+00 
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0.0575 

0.2115E+00 
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0.0675 
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0 . 7528E+00 
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0.0700 
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-0. 2853E+  00 

0 . 9281E+00 

2.175 

0.0725 

0. 1002E+01 

0. 2734E-03 

0 . 1 001E+01 

2.250 

0.0750 

0. 9960E+00 

0.2938E+0Q 

0 • 9538E+00 

2.325 

0.0775 

0. 9248E+00 

0 . 54  03E  +  00 

0 • 7951E+00 

2.400 

0.0800 

0. 7987E+00 

0. 6977E+00 

0. 5584E+00 

2.475 

0.0825 

0. 6365E+00 

0.741 9E+00 

0 . 2933E+00 
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0.0850 

0.4615E+00 

0.6772E+00 

0 . 5312E-01 
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0.0875 
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0.0900 
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0.0925 

0. 6712E-01 
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0.0950 
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0.811 2E- 01 
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0.1050 
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0.1 645E+00 

3.225 

0.1075 
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0 . 2206E+00 

0 . 1226E+00 

3  .300 
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0 . 5233E-01 
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0.1125 

0.5107E-01 

0 . 2252E  +  00 
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0.1150 
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0. 1703E+00 

-0 . 6794E-01 
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0.1175 

0.1642E-01 
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3.600 
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0.5425E-02 

-0. 2201E-02 

3. 750 
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0 . 5098E-01 
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0 • 5124E-02 

4  .200 

0.1400 
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-0 . 3025E-01 

4.275 
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-0. 4311E-01 
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0.1450 

FIG. 6.  FOUR IER  TRANSFORM  OUTPUT  OF  FIG. 5.  DATA 
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FOURIER  TRANSFORM  PLOT  OF  FIG. 5.  DATA 


In  like 


Using  (4.7),  F  (k)  =  (N/2)  A  cos  0  at  k  =  k'. 

fashion  F  (k)  =  (N/2) A  sin  0  at  k  =  k'.  So 
s 

0  =  arctan(F  /F  )  (4.8a) 

s  c 

3* 

and  A  =  ^  (F2  +  F2)  2  .  (4.8b) 

N  c  s 

The  amplitude  and  phase  of  the  modulation  can  then  be 
determined  from  the  computer  printout  values  of  SUMSIN, 
SUMCOS ,  and  SUMSQD  at  a  maximum  in  the  graph  of  the 
transform. 

Several  factors  contribute  errors  to  an  amplitude 
and  phase  determined  in  this  manner.  A  glance  at 
Fig.  6  shows  that  the  phase  of  the  beat  is  0°  at 
2.25  GHz  but  is  -17°  and  +18°  at  adjacent  printed 
values  corresponding  to  Am  =  ? 1/10;  this  is  for  synthe¬ 
tic,  noise  free  data.  The  analysis  is  phase  sensitive. 
This  coupled  with  an  inherent  uncertainty  in  the  foil 
position  makes  absolute  phase  determination  with  the 
transform  method  uncertain  by  about  50°. 

Another  factor  which  contributes  difficulties 
affects  the  amplitude  for  low  frequencies.  By  low  is 
meant  perhaps  2  to  7  complete  beats  in  a  data  sample. 
The  background  fit  can  be  sufficiently  distorted  by 
the  presence  of  beats  to  reduce  the  amplitude  by  as 
much  as  20%.  Data  was  synthesized  (noise  free)  with 
a  single  beat  of  known  amplitude  and  frequencies 


' 
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corresponding  to  2,  3,  4,  5  and  7  complete  oscillations 
over  the  data  sample.  The  data  was  then  transformed 
twice;  first  with  the  known  background  (same  exponen¬ 
tial  as  before)  then  after  being  HOMER  fitted.  This 
was  done  separately  for  a  cosine  and  sine  beat.  The 
results  for  the  HOMER  fit  are  shown  in  Table  1.  In 
all  cases  the  amplitude  error  using  the  true  background 
function  was  less  than  0.2%.  The  last  two  entries  in 
the  table  are  for  a  data  sample  with  two  beats,  one  a 
sine  the  other  a  cosine.  In  all  cases  the  cosine  beat 
had  a  much  larger  amplitude  error  than  the  sine  for 
equal  frequencies.  The  error  also  decreases  with  highe 
frequency  beats  to  a  minimum  of  about  2%. 


Table  1 

Beat  Amplitudes  After  HOMER  Fit 
(True  amplitude  is  5%) 


Frequency 

(GHz) 

Number  of 
Oscillations 

Beat  Amps,  ii 
Cosine 

i  Percent 
Sine 

1.5 

2 

3.96 

4.77 

2.25 

3 

4.17 

4.95 

3.0 

4 

4.45 

4.91 

3.75 

5 

4.59 

4.91 

5.25 

7 

4.75 

4.90 

3.0/5.25 

4/7 

4.30 

4.60 

5.25/3.0 

7/4 

4.76 

4.62 
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Another  source  of  significant  error  to  the  beat 
amplitude  is  due  to  averaging  over  a  finite  length  of 
beam.  This  discussion  will  begin  with  a  single 
modulation  superposed  on  a  background  of  one  exponen¬ 
tial.  It  will  then  be  expanded  to  include  a  background 
of  any  number  of  exponentials  and  the  correction  factor 
so  obtained  can  be  applied  to  any  number  of  modulations. 

Let  the  assumed  intensity  be  of  the  form 

y(x)  =  Bk(x)  A  cos(ax-0)  +  Bk(x)  (4.9) 


with  Bk(x)  =  B  exp (-bx) .  It  is  somewhat  more  convenient 
to  use  the  variable  x  rather  than  the  time,  t,  since 
foil  excitation.  Then 

X  =  vBt 


2ttv 


b  =  1/ (vbt) 


with  v  the  frequency  of  modulation  and  j  the  lifetime 
of  the  excited  state.  Then  the  integrated  intensity  is 


I(x) 


x+A 

» 

y (x' )  dx 1 

x-A 


(4.10) 


with  2A  being  the  observed  beam  length  which  is  dictated 
by  the  geometry  of  the  system.  First  consider  the  inte¬ 
gral  of  the  background: 


■ 
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x+A 

|  Bk(x')  dx'  =  B' e_bx  e  Bk' (x) 
x- A 

where 

_  2B  sinh(bA) 

B - E -  ' 

The  integral  of  the  first  term  of  (4.9)  is: 

x+A 

» 

Bk(x')  A  cos (ax' -0)  dx'  = 

x-A 


2b  e  A 

- ~ - * — {cos(aA)  sinh(bA)[b  cos  (ax-0 )- a  sin  (ax-0 ) ] 

(a^+b^) 

+  sin(aA)  cosh(bA)[b  sin(ax-0)+a  cos (ax-0)]}. 

This  can  be  simplified  and  expressed  in  terms  of  B' 
as 


— hv 

I  (x )  —  Bk'  (x)  =  B'  e  C (a,b ,  A)  A  cos  (ax-0  1  )  , 


(4.11) 


with 


i_  [cos2  (aA)  +  sin2  (aA)  coth2  (bA)  ] 

C  (a ,b  ,  A) - 5 - o— e - 

[1+  (aA)  V(bA)^]  2 


(4.12) 


and  0'  =  0  +  6  +  e.  These  phase  factors  are 
tan  6  =  b/a 

(4.13) 

tan  e  =  -cot(aA)  tanh(bA)  . 
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I (x)  represents  the  data  while  Bk ' (x)  is  the 
background  fit.  As  explained  previously  the  background 
is  subtracted  from  the  data  and  this  difference  is 
divided  by  the  background.  C(a,b,A)  is  then  the  correc¬ 
tion  factor  to  the  amplitude.  For  both  aA  and  bA  <<  1, 

C  =  1  as  expected.  The  condition  that  aA  be  small  is 
that  the  averaging  length  be  a  small  fraction  of  the 
beat  wavelength.  For  bA  to  be  small  the  exponential 
background  must  not  decay  too  rapidly.  Ideally  both 
of  these  conditions  should  be  met  but  in  practice, 
especially  with  non-focused  optics,  this  is  generally 
not  the  case.  Before  considering  some  typical  values 
let  the  background  be  a  sum  of  exponentials: 

Bk  (x)  =  l  Br  exp(-bnx)  . 
n 

Then 

x+A 

Bk(x')  =  £  b'  exp  (-b  x)  =  Bk '  (x) 
x-A  n 

with 

2B  sinh(b  A) 

B’  =  — - k - —  • 

n  b 

n 

For  this  case 

-b  x 

I (x) -  Bk' (x)  =  l  B^  e  n  Cn(a,bn,A) A  cos(ax-0') ,  (4.14) 

n 


with  CR  defined  as  before  with  b  replaced  by  bn»  Here 
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the  added  condition  that  0=0'  has  been  used.  It 
will  be  shown  subsequently  that  this  holds  approxi¬ 
mately  under  "usual"  experimental  conditions.  (4.14) 
shows  that  normalization  hopelessly  muddles  the  picture 
unless  all  the  Cn  are  approximately  equal.  This  is  in 
fact  the  case  for  aA  <  tt/2  but  differences  rapidly 
become  too  large  to  ignore  as  aA  approaches  tt. 

If  in  addition  there  is  more  than  one  modulation 
the  expression  for  the  residuals  becomes 

-b  x 

I(x)-Bk'(x)=  y  b'  e  n  C  (a  .b  ,A)  A  cos (a  x-0') 

L  n  nm  m'  n'  m  m  m 

n,m 


(4.15) 


and  a  correction  factor  exists  for  each  A  . 

m 

Since  2 A  is  the  averaging  length,  aA  =  tt  corres¬ 
ponds  to  a  length  of  one  wavelength  and  the  condition 
that  aA  <  tt/2  is  almost  mandatory  for  a  meaningful 
correction  to  the  beat  amplitude.  For  this  experiment 
A  =  0.6  mm  which  gives  aA  =  1.1  for  hydrogen  and 
aA  =  1.4  (typically)  for  boron  V.  For  hydrogen  which 
has  a  moderate  single  exponential  background,  bA =  .05 
while  for  boron  V  with  a  steep  double  exponential 
b-^A  =  1.0  and  b2A  =  0.25.  With  these  numbers  it  is 
possible  to  show  the  behavior,  under  somewhat  extreme 
conditions,  of  C(a,b,A)  and  0'. 


. 
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The  phase  shift  introduced  by  averaging  over  a 
finite  length  of  beam  is  0-0'  =  6  +  e.  From  (4.13) 

lim  (6+e)  =  0 
aA+0 

regardless  of  bA.  For  aA  =  tt  this  phase  shift  is  about 
tt/2  for  any  value  of  bA.  Table  2  illustrates  some  inter¬ 
mediate  values.  That  0'  ~  0  for  aA  <  tt/2  and  bA  <  1  is 
seen  to  hold. 


Table  2 

Beam  Averaging  Phase  Shift 


aA 

CD 

1 

CD 

bA  =  .05 

bA  =  .25 

bA  =  1.0 

0 

0 

0 

0 

it/  4 

O 

00 

• 

O 

4° 

15° 

tt/2 

2° 

9° 

32° 

it 

91° 

95° 

108° 

Fig.  8  is  a  set  of  graphs  of  C(a,b,A)  vs  aA  for 
three  different  values  of  bA  (the  same  ones  as  in 
Table  2)  .  For  aA  <  tt/2  C  is  insensitive  to  bA  justi¬ 
fying  its  use  as  a  single  correction  factor  even  for  a 
multi-exponential  background.  For  aA  >  tt/2  this  is  not 
the  case  and  a  straightforward  correction  to  the  beat 
amplitude  is  not  possible. 
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Fig.  8  Beat  Amplitude  Correction  Factor 
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The  analysis  considered  so  far  is  still  deficient 
for  determining  absolute  values  of  beat  amplitude.  The 
integrated  intensity  as  given  by  (4.10)  should  be 
replaced  by 


I  (x) 


x+A 

R  (x  ' )  y  (x ' )  dx 1 

o 

x-A 


with  R  the  proper  response  function  of  the  experimental 
apparatus.  However  frequently  only  relative  amplitudes 
of  several  beats  appearing  in  the  Fourier  transform  are 
desired.  For  this  (4.15)  with  the  modified  version  of 
(4.12)  supplies  the  necessary  caution  that  the  relative 
amplitudes  appearing  in  the  transform  are  not  the 
relative  amplitudes  of  the  beats  as  well  as  an  expression, 
valid  under  broad  but  restricted  conditions,  to  correct 
for  this. 


4 . 3  Least-Squares  Evaluation  of  Amplitude  and  Phase 

A  second  method  to  determine  the  amplitude  and 
phase  of  a  modulation  superposed  on  a  multi-exponential 
background  was  developed  using  a  straightforward 
linear  least-squares  fit.  For  linearity  the  frequency 
is  held  constant  and  input  at  the  value  given  by  the 
Fourier  transform,  but  the  program  is  inexpensive 
enough  to  scan  a  selection  of  frequencies  and  produce 
a  plot  much  like  the  transform  method.  The  advantage 
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over  the  transform  method  is  that  uncertainties  in  both 
amplitude  and  phase  can  be  evaluated  and  the  reduced 
chi-squared  with  and  without  a  modulated  background 
can  be  compared.  At  present  this  program  is  limited 
to  a  single  frequency  modulation  but  expansion  to 
include  up  to  five  frequencies  (10  parameter  fit)  should 
not  be  difficult  (and  should  be  completed  prior  to  the 
next  experiment) . 

The  method  of  least-squares  is  to  minimize  chi- 
squared,  the  weighted  sum  of  the  squares  of  the  devia¬ 
tion  between  the  data  and  the  fitting  function.  That 
is. 


N 


I 

i=l 


y(V]2 


with  the  datum  at  x^ ,  y(x^)  the  fitting  function, 
and  o ^  the  uncertainty  in  the  datum  which  for  photon 
counting  experiments  is  taken  to  be  the  square  root  of 
the  datum,  assuming  Poisson  statistics.  The  fitting 
function  is  given  by 


y (x . )  =  Bk  (x . )  +  Bk (x . )  [A, cos (ax . ) +  A~sin (ax . ) ]  . 

J  1  l  i  1  i  2  l 


Here  Bk(x^)  is  the  background  fit  as  determined  by  the 
HOMER  routine,  a  is  proportional  to  the  assumed  known 
frequency  and  A^  and  A£  are  parameters  to  be  determined. 
For  these  the  amplitude  and  phase  relative  to  a  cosine 
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modulation  are  easily  determined  from  elementary 
relations . 

2 

Setting  the  partial  derivatives  of  x  with 
respect  to  A^  and  A ^  equal  to  zero  leads  to  the 
linear  equations: 


f  . 
3 


2 

I 

k=l 


a  ..  A, 
Dk  k 


with 


N 


f  . 
3 


=  I 


i=l 


a . 

i 


(yiBk(xi)-  Bk(xi)2) 


cos (X  .  )  ^ 
1 

^■sin  (x^ ) 


and  j  is  the  column  index  while 

2  ,  x  . 
cos  (x^)  cos.sm 

2 

cos • sin  sin  (x^ ) 

Written  this  way  it  is  clear  how  to  expand  the  technique 

to  more  than  2  parameters.  Note  that  it  is  the  data 

itself  which  is  fit  and  not  the  normalized  residuals. 

In  this  way  a  ^  is  the  error  matrix  and  the  uncorrela- 

2  -1 

ted  uncertainties  in  A.  are  just  aa  =  a. .  (Be  69). 

9  A  j  3  3 

A^  and  A^  are  obtained  in  the  usual  way  and 
these  are  readily  converted  to  a  cosine  amplitude  and 
phase.  The  uncertainties  in  the  parameters  A  and  0 


N 


a 


jk 


=  l 


Bk (x. ) 
1 


i=l 


with 


. 
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A  cos  (ax-9)  =  A-^cos  (ax)  +  A2sin(ax) 


A 


2 


A 


2 


+  A 


2 


1 


2 


are  related  to  the  uncertainties  in  A^  and  A by 


and  similarly  for  9 


Fig.  9  shows  an  abbreviated  output  of  this 
program  for  the  same  synthetic  data  as  Fig.  5. 

The  most  obvious  difference  between  them  is  that  the 
graph  of  the  least-squares  output  does  not  go  to  zero 
at  1.5  and  3.0  GHZ  since  there  is  no  orthogonality 
condition  for  this  method.  The  beat  frequency  is 
decided  by  a  minimum  in  the  reduced  chi-squared  rather 
than  the  position  of  maximum  amplitude  since  these  two 
criteria  do  not  necessarily  agree.  Reduced  chi-squared 
is  simply  chi-squared  divided  by  the  number  of  degrees 
of  freedom  for  the  fit.  This  number  is  the  number  of 
data  points  less  the  number  of  fitting  parameters  which 
should  include  the  number  of  HOMER  parameters  even 
though  these  are  not  being  varied.  How  close  the 
reduced  chi-squared  is  to  one  is  an  indication  of  the 
goodness-of-f it .  Using  this  criterion  one  can  obtain 
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the  "best"  values  for  frequency,  amplitude  and  phase 
with  uncertainties  in  each.  It  must  be  mentioned 
that  if  more  than  one  frequency  were  present  on  the 
data,  even  widely  separated  frequencies,  the  answers 
obtained  with  this  method  would  be  incorrect  since 
the  reduced  chi-squared  would  never  be  close  to  one. 

A  least-squares  fit  with  as  many  frequencies  as  are 
present  on  the  transform  must  initially  be  used.  The 
complementary  nature  of  these  two  techniques  thus 
becomes  apparent.  While  the  transform  method  gives 
an  indication  of  the  amplitude  and  phase  for  several 
frequencies  it  offers  no  indication  of  the  reliability 
of  the  results.  Least-squares  gives  a  clear  criterion 
for  goodness-of-f it  but  becomes  increasingly  complicated 
as  more  than  one  frequency  are  fit. 
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CHAPTER  V 


RESULTS  AND  CONCLUSION 


Before  the  analysis  of  any  beam-foil  data  involv¬ 
ing  time-of-f light  information,  the  beam  velocity  must 
be  determined.  Two  methods  are  usually  used.  The 
first  relies  on  the  energy  calibration  of  the  accele¬ 
rator  and  compensates  for  energy  loss  in  the  foil  using 
published  tables  (NS  70).  The  second  is  a  direct 
independent  measurement  of  the  velocity  usually  invok¬ 
ing  the  Doppler  shift  of  a  spectral  line.  This  involves 
a  specially  designed  target  chamber  with  an  arrangement 
of  mirrors  which  has  been  described  elsewhere  (Li  71) . 
The  first  method  was  applied  to  all  elements  examined 
during  the  Spring  1974  experiment  at  Alberta  since 
previous  Doppler  shift  measurements  indicated  an 
accuracy  of  the  beam  velocity  to  within  2%.  The  energy 
loss  for  hydrogen  and  boron  at  the  beam  energies  employed 
was  less  than  1%  and  this  correction  was  never  applied 
to  the  quantum  beat  data  since  it  was  not  justified  by 
the  results.  Highly  accurate  quantum-beat  frequency 
determination  can  only  be  obtained  by  comparing  two 
observed  frequencies,  one  of  which  is  assumed  accurately 
known  from  other  theoretical  or  experimental  evaluation. 
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5.1  Electric  Field  Calibration 

The  target  chamber  at  Argonne  National  Laboratory 
was  designed  specifically  to  yield  a  uniform,  longitu¬ 
dinal  electric  field  along  the  beam  axis.  Calibration 
was  accomplished  by  measuring  the  frequency  of  quantum 
beats  observed  in  hydrogen  Lyman-a  emission  (ls-2p  at 
1216  A).  H*  was  accelerated  to  1.0  MeV  and  a  double 
foil  was  used,  the  first  to  dissociate  the  molecule 
into  three  protons.  Fig.  10  is  a  semi-log  plot  of  the 
data  at  -1,000  volts  and  Fig.  11  is  its  Fourier  trans¬ 
form.  Measurements  were  taken  at  three  different  plate 
voltages  and  these  with  the  observed  frequencies  and 
their  uncertainties  are  given  in  Table  3.  The  uncer¬ 
tainties  are  solely  the  statistical  estimates  given  by 
the  least-squares  routine.  The  2%  uncertainty  in  the 
beam  velocity  is  greater  than  any  of  these.  A  calibra¬ 
tion  curve  for  Ly-a  was  calculated  as  mentioned  in 
section  2.6.  The  field  strengths  for  these  three  points 
are  shown  in  Table  3.  The  first  two  points  agree  that 
the  electric  field  is  6%  higher  than  the  plate  voltage 
divided  by  the  2  cm  separation.  The  agreement  indi¬ 
cates  a  systematic  error  which  can  be  fully  accounted 
for  by  the  beam  velocity  uncertainty  and  the  power 
supply  voltage  tolerance. 
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Fig. 10  Actual  Hydrogen  Lyman-a  Data 
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Table  3 

Quantum  Beat  Frequencies  for  Lyman-g 


Plate  Voltage 
(Grounded  Foil) 

Frequency 

(GHz) 

Uncertainty 

Field  Strength 
(volts/cm) 

-1,000 

2.44 

±.01 

530  ±  3 

-2,000 

4.53 

±.05 

1060  ± 14 

-4,000 

(8.0) 

(±.l) 

(1970  ±  30) 

The  figures  for  -4,000  volts  appear  in  paren¬ 
thesis  since  the  beat  amplitude  was  so  small  that  it 
hardly  appeared  on  the  Fourier  transform.  The  least- 
squares  routine  gave  the  tabulated  frequency  with  a 
reduced  chi-squared  of  1.02  compared  to  1.15  without 
any  modulation.  The  small  beat  amplitude  is  a  result 
of  averaging  over  one  beat  length.  At  this  frequency 
aA  =  3.76  and  Fig.  8  claims  that  the  beat  amplitude 
will  be  reduced  over  90%. 

5 . 2  Results  with  Boron 

Boron  III  2066  A 

2  2 

The  2s  Si  -  2p  Pi  3  transition  of  boron  III  at 

,  -55 

2066  A  was  examined  for  zero-field  quantum  beats  at 
0.6  and  1.0  MeV.  and  ^B  were  individually 

analyzed  after  taking  data  samples  with  80  points 
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separated  by  1.0  mm.  No  polarizer  was  used  and  the 
line  was  sufficiently  intense  to  collect  good  statis¬ 
tics.  The  decay  was  a  single  exponential  with  a 
lifetime  of  5.5  ns  which  at  the  lower  beam  energy 
corresponds  to  18  channels.  has  a  nuclear  spin, 

1=3  while  "^B  has  I  =  3/2.  The  calculated  hyperfine 
structure  frequencies  fall  well  within  experimental 
resolution  but  no  beats  were  consistently  detected  at 
both  beam  energies.  Also  no  observed  beats  had  rela¬ 
tive  frequencies  consistent  with  theory. 

It  might  be  premature  to  conclude  that  no  beats 
are  on  the  data.  Since  multiple  frequencies  are 
expected,  the  least-squares  program  could  not  be  used. 
Furthermore  only  F  value  intervals  corresponding  to 
J  =  3/2  were  calculated.  Since  the  fine  structure  is 
not  optically  resolved  (34.1  cm  and  is  small  enough 
for  coherence  to  be  present,  J  =  1/2  should  also  be 
considered.  The  theoretical  frequency  pattern  could 
be  much  more  complicated  than  the  one  we  employed. 

This  data  will  be  among  the  first  analyzed  by  the 
expanded  version  of  the  least-squares  routine  and  if 
the  analysis  clearly  indicates  the  existence  of  beats 
a  detailed  theoretical  investigation  of  the  hyperfine 


structure  will  be  undertaken. 
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Boron  V  749  A 

Data  at  749  A  (n  =  4  to  n  =  3)  for  ^B  V  was 
collected  at  Argonne  National  Laboratory  at  7  different 
plate  voltages  between  -1,000  and  -5,000  volts.  Beam 
energies  of  3.2  or  3.0  MeV  were  used.  Typically  each 
data  sample  represents  six  scans  along  the  beam  with 
43  points  per  scan.  Each  point  was  spatially  separated 
by  0.27  mm  which  represents  a  temporal  separation  of 
about  .04  ns.  Each  scan  was  added  with  the  on-line 
ASI  2100  computer.  Table  4  lists  some  of  the  relevant 
parameters  for  each  data  sample.  The  lowest  and  highest 
plate  voltages  were  taken  last  at  the  lower  beam  energy 
in  an  effort  to  improve  the  statistics  by  using  higher 
beam  current. 


Table  4 


Experimental  Parameters  for  Boron  V  749  A 


Plate  Volt, 
(volts ) 

Beam  Vel. 
(mm/ns ) 

Scale  Factor 
(cm) 

Aver . Length 
(mm) 

Statistics 

(%) 

-1,000* 

7.25 

1.13 

1.2 

2.4 

-2,000 

7.49 

1.10 

2.0 

2.0 

-2,500 

7.49 

1.10 

2.0 

2.0 

-3,000 

7.49 

1.10 

2.0 

2.0 

-3,500 

7.49 

1.10 

2.0 

2.0 

-4,000 

7.49 

1.10 

2.0 

2.0 

-5,000* 

7.25 

1.13 

1.2 

2.2 

* 

Beam  energy  3.0  MeV,  otherwise  3.2  MeV. 
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The  fourth  column  of  Table  4  lists  the  beam 
averaging  length  as  described  in  Chapter  III.  The 
first  and  last  data  sample  used  a  larger  grating 
mask  and  smaller  slit  width  than  did  the  others. 

The  frequency  whose  wavelength  equals  the  averaging 
length  is  3.75  GHz  for  2.00  mm  and  6.04  GHz  for  1.2  mm 
at  the  listed  velocities.  These  represent  the  prac¬ 
tical  upper  frequency  limits  which  are  considerably 

smaller  than  1/(2S  )  t  13.5  GHz.  Counting  time  for 

c 

the  two  data  samples  with  higher  resolution  was  twice 
as  long  as  for  the  others  and  still  fewer  counts  per 
datum  were  collected.  This  is  represented  by  the 
"Statistics"  column  which  is  the  mean  value  over  a 
data  sample  of  the  square  root  of  a  count  divided  by 
the  count. 

A  simplistic  but  reasonable  assertion  might  be 
that  one  cannot  detect  beats  of  amplitude  less  than 
the  statistics  of  the  data  sample.  As  a  generality 
this  is  not  true  since  the  noise  level  of  a  Fourier 
transform  depends  on  how  many  data  points  have  been 
collected.  However,  for  43  points,  the  assertion 
appears  correct  based  on  studies  of  synthesized  data 
which  were  "fuzzed"  as  described  in  the  previous 
Chapter.  Figs.  12  to  14  show  the  Fourier  transforms 
of  synthetic  data  with  three  widely  separated  frequen¬ 
cies.  The  data  was  carefully  generated  to  duplicate 
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FIG. 12 


FOURIER  TRANSFORM  OF  SYNTHETIC  DATA  I 
(FREQS  IN  GHZ:  1.8,3. 6,5. 4;  AMP:1%  NO  FUZZ) 
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FIG. 13.  FOURIER  TRANSFORM  OF  SYNTHETIC  DATA  II 

(FREyS  IN  GHZ;  1.8, 3. 6, 5. 4;  AMP:2%  FUZZ) 
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FIG. 14.  FOURIER  TRANSFORM  OF  SYNTHETIC  DATA  III 

(FREQS  IN  GHZ:  1.8, 3. 6,5.4;  AMP : 1  %  FUZZ) 
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actual  boron  data  with  2%  statistics.  The  first  has  a 


1%  beat  amplitude  with  no  scatter.  The  second  and 
third  have  artificial  scatter  with  beat  amplitudes  of 
2%  and  1%  respectively.  The  2%  case  has  been  strongly 
effected  by  FUZZ  but  the  three  frequencies  are  clearly 
larger  than  the  noise.  For  the  1%  case  the  "true" 
frequencies  cannot  be  distinguished  from  the  background. 

Fig.  15  is  a  graph  of  theoretically  possible 
frequencies  calculated  with  the  Stark-effect  program 
mentioned  in  Chapter  II.  Three  curves  converge  at  zero 
field  to  the  P-D  Lamb  shift  at  1.42  GHz.  Five  converge 
to  the  D-F  Lamb  shift.  No  theoretical  frequency  is 
higher  than  the  upper  curve  and  lower  than  51  GHz,  the 
value  of  the  S-P  Lamb  shift  at  zero  field.  It  was 
hoped  that  Stark  beats  could  be  observed  that  would 
follow  one  or  more  of  the  upper  curves  to  the  P-D  Lamb 
shift.  In  fact,  no  beat  was  observed  with  an  amplitude 
greater  than  the  statistics  value  of  the  data  sample. 

In  all,  33  beats  were  observed  with  amplitudes  between 
1%  and  2.4%,  and  all  fall  randomly  between  3  and  13.5 
GHz.  The  conclusion  is  that  these  beats  are  noise  and 
Stark  beats  in  Boron  V  were  not  observed. 

5 . 3  "Eck  Beats" 

In  addition  to  the  three  data  samples  of  hydro¬ 
gen  taken  for  field  calibration,  a  fourth  set  of  data 
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was  taken  at  +1,000  volts.  The  minus  and  plus  1  KV 

data  allow  the  relative  cross  section  for  S-P 

coherence,  a  ,  to  be  evaluated,  as  in  section  2.6. 
sp 

The  calculation  will  be  outlined  and  the  results 
presented  more  as  an  indication  of  current  quantum 
beat  research  rather  than  the  validity  of  the  results. 

The  beat  amplitudes  were  corrected  for  finite 
beam  acceptance  length  and  the  beats  were  added  and 
subtracted  yielding: 

I+  +  I_  =  .100  cos  ( oot  -  80°)  (5.1) 


I  -  I_  =  .056  cos  (cot  +  36°) 


(5.2) 


To  proceed  with  the  calculation  the  values  of  a  ,  a 

s  Po 

and  a  are  required  and  these  were  taken  from  GAM  74. 
Pi 

Since  only  the  relative  values  are  known,  the  sum 

a  +  a  +  2a  has  been  normalized.  a  =  0.56, 
s  p0  Pi  s' 

a  =  0.14  and  a  =  0.15.  If  equation  (5.1)  is 
Po  Pi 

compared  to  (2.11),  the  modulation  can  be  made  a 
negative  cosine  by  adding  -100°-  to  the  phase  of  (5.1). 
This  is  equivalent  to  shifting  the  first  point  of  the 
sum  back  to  the  foil  position.  The  constant  term,  V, 
may  then  be  evaluated.  By  adding  the  same  phase  factor 
to  (5.2)  and  using  V,  (2.12)  may  be  used  to  evaluate 


a 


sp 
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Table  5  lists  the  results  for  this  experiment 
at  333  KeV  and  those  of  Gaupp  et  al.  (GAM  74)  at  the 
other  beam  energies.  The  discrepancy  motivated  a 
repeat  of  this  part  of  the  experiment  by  the  group  at 
Argonne  National  Laboratory  without  the  author's  par¬ 
ticipation.  Attention  was  focused  on  duplicating  the 
starting  point  of  the  foil  by  taking  plus  and  minus 
data  samples  alternately  and  on  obtaining  better 
statistics.  The  analysis  of  this  new  data  is  not  yet 
complete . 


Table  5 

Cross  Section  for  S-P  Coherence 


H+  Energy 

Re  o 

sp 

Im  o 

sp 

1  G 

1  sp1 

110  KeV 

-0.03 

0.06 

0.07 

160  ” 

-0.09 

0.17 

0.19 

210  " 

-0.10 

0.19 

0.22 

333  " 

+0.14 

0.12 

0.19 

5 . 4  Conclusion 

To  observe  quantum  beats  many  necessary  condi¬ 
tions  must  be  satisfied.  Chapter  II  discussed 
coherence  and  alignment.  The  former  is  guaranteed 
by  the  thin  foil  but  the  latter,  while  known  to  be 
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a  property  of  foil  excitation,  is  not  guaranteed  to  be 
of  sufficient  magnitude  to  produce  observable  beats 
for  any  sample.  Furthermore  the  mechanism  for  align¬ 
ment  is  not  yet  understood  so  it  is  not  possible  to 
forecast  the  degree  of  alignment  theoretically.  Good 
spatial  resolution  is  an  obvious  requirement  for  obser¬ 
vation  of  high  frequencies.  When  experiments  are 
performed  with  equipment  pushed  to  its  design  limit, 
as  was  the  case  with  hydrogenic  B  V  line  in  this 
project,  it  is  easy  not  to  satisfy  the  requirement 
and  be  unaware  of  this  until  after  the  experiment  has 
been  completed.  The  requirement  for  good  statistics 
cannot  be  overemphasized.  Often  the  spectral  line  of 
theoretical  interest  is  not  very  intense  and  a  com¬ 
promise  struck  in  the  laboratory  proves  inadequate  on 
the  computer.  The  effect  of  cascades  must  also  be 
considered.  Whether  or  not  alignment  can  be  transferred 
by  a  cascade  has  yet  to  be  definitely  established.  It 
seems  best  to  choose  a  cascade-free  transition  if  at 
all  possible. 

It  is  too  early  to  conclude  that  the  2066  A 
measurement  was  unsuccessful.  Hyperfine  beats  are 
noted  for  their  small  amplitudes  so  further  analysis 
is  desirable.  There  are  no  quantum  beats  on  our  boron 
V  data.  To  be  properly  performed  the  experiment 
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requires  better  spatial  resolution  and  good  statistics, 
which  might  be  obtained  if  the  4  to  2  or  4  to  1  tran¬ 
sition  is  observed  with  a  grazing  incidence  spectro¬ 
meter,  which  we  do  not  have.  The  hydrogen  quantum 
beats  were  clearly  observed  and  the  "Eck  beat"  calcu¬ 
lation  which  has  been  presented  is  both  current  and 
controversial . 

This  work  has  been  an  attempt  to  review  the 
entire  field  of  quantum  beat  research  and  to  present 
a  description  of  early  experiments  and  analysis  per¬ 
formed  at  the  University  of  Alberta.  It  is  hoped  that 
future  work  here  will  benefit  from  this  endeavor. 
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FOURIER  TRANSFORM  USING  HOMER  FIT 

DIMENSION  OSCN  (400)  ,  BEAT  (400)  ,  BKGD  (400)  ,  LINE  (71) 
INTEGER  DO I, STAR, BLANK, SIG (400) 

INTEGER  WORDS  (72) 

DATA  DOT/1 H ./, SIAR/1 H*/ ,  BLANK/1  H  / 

50  READ  (5,201) N , MIN , MAX , SCALE , STEP, FINC R , SIZE, ENLO 
IF(N.LT.I)  GO  TO  60 
READ  (5,205)  WORDS 
READ  (4 , 202 )  (SIG  (I)  ,  1  =  1 ,  N) 

READ  (3,210)  IA , A0 , A2  ,TB, A4 , TC , VEL 
WRITE (6,100)  WORDS 
WRITE  (6, 109) SCALE 

EVALUATE  BACKGROUND  PARAMETERS 
USING  HOMER  FIT 

IF  (TA)  16,16,17 

16  A 1 =0 . 

GOTO  20 

17  A1=-STEP/(VEL*TA) 

IF(TB)  20,20,2  5 

20  A3  =0 . 

GOTO  30 

25  A3=-STEP/(VEL*TB) 

IF  (TC)  30,30,35 
30  A5=0 . 

GOTO  40 

35  A5=-STEP/(VEL*IC) 

CALCULATE  BEAT  AMPLITUDE 

40  BTSUM=0 . 0 
DO  6  1=1, N 
QI=(I-D 

BKGD  (I) =A0*EXP  (A1*£I) +A2*EXP (A3*QI) 

*  +  A4*EXP  (A5*QI) 

OSCN  (I)  =SIG  (I)  -BKGD  (I) 

IF  (BKGD  (I))  5,5,7 

5  BEAT  (I) =OSCN  (I) 

GOTO  6 

7  BEAT  (I) =OSCN  (I) /BKGD  (I) 

6  BTSUM=BTSUM  +  BEAI (I)  *BEAT (I) 

APRX=8 . 0*BTSUM*N/9. 8696 
WRITE  (6 , 110)  APEX 

WRITE (6,104) 

WRITE  (6, 102)  (SIG  (I)  ,  BKGD  (I)  ,  OSCN  (I)  ,  BEAT  (I)  ,1=1,  N) 
CALCULATE  FOURIER  TRANSFORM 


; 
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WHITE  (6,101) 

DO  15  1=1,71 
15  LINE  (I)  =BLANK 
DO  10  I  =  MIN ,  MAX 
LINE  (  1)=DOI 
QS=I-1 

EN=ENLO+FINCR*QS 
SUMSIN=0.0 
SUMCOS=0 .0 
P=EN*SCALE 
DO  11  J=1,N 
QT=J-1 

ARG=2.*3.14159*P*QT 

SU M SIN  =  SUM SIN*  BEAT (J) *SIN  (AEG) 

SUMCOS=SUMCOS+BEAI  (J) *COS  (ARG) 

11  CONTINUE 

SUMSQD=SUMSIN*SUMSIN+SUMCOS*SUMCOS 
FREQ=EN*30 . 0 

PLOT  FOURIER  TRANSFORM 

IY  =1 . +SUMSQD*SiZE/APRX 
LINE  (IY) =STAR 

WRITE (6, 107) SUMSQD,SUMSIN, SUMCOS, FREQ, EN, LINE 
LINE (IY) =BLANK 
10  CONTINUE 
GO  TO  50 
60  WRITE (6 , 108) 

STOP 

100  FORMAT (1H1 ,72A1) 

101  FORMAT  ( » 1 ■ ,  4X  , 'SUMSQD • ,6X, ' SUM SI N* ,6X,*SUMCOS'  ,4X, 
*'FREQU. • ,2X, 'ENERGY' ,30X, ’AMPLITUDE  OF  FOURIER', 

*•  TRANSFORM' ,/39X,'  (GHZ)  ' ,3X, '  (CM-1)  ' , 

*7X,  *  0 . 1 . 2 . 3 . •  , 

*'4 . 5 . 6 . 7') 

102  FORMAT (110, 5X,3E15. 6) 

104  FORMAT (1  HO, 4X, 'SIG (I)  ' ,8X,'BKGD  (I)  ' ,8X,'OSCN(I)  * 

*  ,  8X,  'BEAT  (I)  '  ) 

107  FORMAT ( 1 H  , 3E1 2. 4 ,2F8 . 4 , 7X, 71 A1 ) 

108  FORMAT ( 1 Hi ,50X, ' THE  END') 

109  FORMAT (1H-, 14HSCALE  FACTOR  =,F6.3, 

*21 H/ (CHANNEL* WAVENUMBER)  ) 

110  FORMAT (1H0, 'MAX  SUMSQD  =  *,E10.4) 

201  FORMAT (313, IX, 5F10. 5) 

202  FORMAT(/,/,  (5  (3X,I6,5X) ) ) 

205  FORMAT (72A1) 

210  FORMAT  (/ ,  1 9X,F9 . 4 , 2F 1 1 . 2 , FI  0 . 4 , FI 0 . 2 , FI  0 . 4 ,/, /, 
*73X,F5. 3) 

END 


' 


..  .«•»%  'V 
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C 

C 


BEAT  AMPLITUDE  EVALUATION  WITH  PLOT 


INTEGER  SIG  (256) , WORDS  (72)  ,DOT, STAR , BLANK ,AMPL (2  56) 
* , FLGPLT 

DIMENSION  LINE  (100)  ,  FREQ  (256)  , BKGD  (256) 

DATA  DOT/1 H ./, SIAR/1 H*/,BLANK/1H  / 

PI=3. 14159 

15  READ  (5  ,  130)  MOLDS 

READ (5,120)  N , MAX, SCALE, STEP, FI  NCR, FREQ1 ,FLGPLT 

IF(N.EQ.O)  STOP 

WRITE  (6,200)  WORDS 

WRITE (6 , 202)  SCALE 

READ  (3, 100)  IDF, NHP 

READ  (3,105)  TA,A0,A2  ,TB , A4 , TC , VEL 

READ  (4,110)  (SIG(I)  ,1  =  1, N) 

A1=-STEP/(VEL*IA) 

IF(TB)  20,20,25 
20  A3=0. 

GOTO  30 

25  A3  =- STEP/ (VEL* IB) 

IF  (TC)  30,30,35 
30  A5=0 . 

GOTO  40 

35  A5=-STEP/(VEL*IC) 


C 

c 

c 


EVALUATE  BKGD,  STATS  AND  CHISQO 


40  CHISQ0=0 . 

STSUM=0 . 

DO  45  1=1, N 

BKGD  (I) =A0*EXP  (A1*  (1-1) ) +A2*EXP (A3* (1-1)  )  + 
*A4*EXP  (A5*  (1-1)  ) 

STSUM=SISUM+SQRI (FLOAT  (SIG (I) ) ) /SIG (I) 

45  CHISQO=CHISQG ♦  (  (SIG  (I)  -BKGD  (I)  )  **2)  /SIG  (I) 
STATS=1 00* STS  UM/N 
NF  P=NHP-IDF 
RFP=FLOAI (N)  /NMP*NFP 
CHISQ0=CHISQ0/ (N-RFP) 

WRITE (6 , 205)  CHISQO,STATS 
WRITE  (6,210) 

WRITE (6 , 220) 

DO  70  L= 1 , MAX 

FREQ  (L) =FREQ1+  (L- 1 ) *FINCR 

Y1=2*PI*FREQ (L) *SCALE/29. 979 


C 

C 

C 

C 


2  PARAMETER  LINEAR  LEAST  SQUARES  FIT 
TO  AMPLITUDE  AND  PHASE 


B1 1  =0 
B 1 2=0 
B2  2=0 
R1  =0 . 


’  J 


•  •  • 


. 


. 
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R2  =0 . 

CHISQ=0. 

ENTER  SUMMATION  LOOP 
DO  50  1=1,  N 

B1=C0S  (XI*  (1-1) ) *BKGD  (I) 

B2  =  SIN  (Y1*  (1-1)  )  *BKGD  (I) 

B11  =  B1H-B1*B1/SIG  (I) 

B12=B12+B1*B2/SIG (I) 

B22=B22+B2*B2/SIG  (I) 

R1=RH-B1-BKGD(I)*B1/SIG  (I) 

R2=R2+B2-BKGD  (I) *B2/SIG (I) 

50  CONTINUE 

EVALUATE  AMPLITUDES  AND  UNCERTAINTIES 

DET=B11*B22-B12*B12 
AMP1  = (R1*B22-R2*B12) /DET 
AM  P2= (B11*R2-R1*B12) /DET 
UAMP1=B22/DET 
UAMP2=B 1 1/DET 
UAMP12=-2*BU/DET 

CHANGE  VARIABLES  TO  COSINE 
AMPLITUDE  AND  PHASE 

AMP1 2=AMP1 *AMP 1 
AMP22= AMP2*AMP2 
AMP=SQRT  (AMP1 2+AMP22 ) 

PHS=ATAN2  (AMP2 / AMP1 ) 

UAMP= (UAMP1*AMP12*UAMP2*AMP22 
*+UAMP12*AMP1*AMP2) / (AMP12+AMP22) 

UPHS= (UAMP1*AMP22+UAMP2*AMP1 2-UAMP12 
+  *AMP1*AMP2)  /(AMP12+AMP22) **2 
UAMP=SQRI (UAMP) *100. 

UPHS=SQRT(UPHS) *360./  (2*PI) 

EVALUATE  CHISQ 

DO  60  1=1, N 

FCTN  =  BKGD  (I) * (1+AMP1 *COS  (Y1*  (1-1) ) +AMP2* 

*  SIN (Y1 *  (1-1) ) ) 

60  CHISQ=CHISQ+ (  (SIG  (I) -FCTN) **2)/SIG  (I) 

AMPP=AMP*1 00. 

PHS=PHS*360 ./  (2*PI) 

CHISQR=CHISQ/  (N-RFP-2) 

AM  PL  (L)  =ABS  (1  000*AMP) 

WRITE  (6  , 230)  FBEy (L)  , AMPP , UAMP, PHS , UPHS , CHISQ ,CHISQR 
70  CONTINUE 

PLOT  RESULTS 

IF  (FLGPLT.EQ. 0)  GOTO  15 


.  '  '  • 


DO  80  1=1,100 
80  LINE (I) =BLANK 
WRITE  (6,250) 

WRITE  (6,260) 

DO  90  1=1, MAX 
LINE  (1)  =DOT 
LINE  (1+AMPL  (I)  )  =STAR 
WRITE  (6,270)  FREQ  (I)  , LINE 
90  LINE  (1+AMPD  (I)  )=BLANK 
GOTO  15 

100  FORMAT  (33X, 13, 31X, 13) 

105  FORMAT  (1 9X, F9. 4, 2F 11. 2, F10. 4, F10. 2, F10. 4, 

*/,  /,  73X,  F5.  3) 

110  FORMAT  (/,/,  (5  (3X,I6,5X)  )  ) 

120  FORMAT (213, 4F10. 5, II) 

130  FORMAT  (7 2A1 ) 

200  F0RMAT(1H1,6X,72A1) 

202  FORMAT (1H0, 6X, 14HSCALE  FACTOR  =,F6.3, 

*21 H/  (CHANNEL*WAVENUMBER)  ) 

205  FORMAT ( 1  HO , 6X , 25HRED0CED  CHISQD  W/O  BEAT  =,F7.3, 

*1 5X, 7HSTAIS  =,F4. 1 ,1 H%) 

210  FORMAT (1H0, 1GX,9HFREQUENCY,7X,9HAMPLITUDE,6X , 

*1 1HUNCERTAINIY,8X,5HPHASE,8X,11H UNCERTAINTY, 

*  1 3X, 1 1 HCHI- SQUARED) 

220  FORMAT (13X,5H  (GHZ) , 9X,9H  (PERCENT) , 9X , 5H  (AMP)  , 9X , 

*9H  (DEGREES)  ,9X,5H  (PHS)  ,25X,9H  (REDUCED)  ,/) 

230  FORMAT (IX, FI 7. 3, FI  6. 3, FI  6. 3, FI  6. 1, FI  6. 2, FI 6. 2  ,F16. 3) 
250  FORMAT (1  Hi , 6X , 9HFREQUENC Y ,42X , 

♦25HBEAT  AMPLITUDE  IN  PERCENT) 

260  FORMAT  (10X,3HGHZ,8X,  20H0 . 1 . . 

*4  0H2 . 3 . 4 . 5 . . 

*42H6 . 7 . 8 . 9 . 10) 

270  FORMAT (1X,F13. 3,7X, 100A1) 

END 


' 
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C  FUZZ  ROUTINE 

IMPLICIT  LOGICAL  (L) , REAL*8  (D) 

INTEGER  STEP  (400) 

C  THIS  SUBROUTINE  RANDOMIZES  DATA  IN  FKS , 

C  FKS  BEING  THE  MEANS 

DATA  GAUSSL  /37.Q/ 

C  FOLLOWING  DEFINES  VARIABLES  FOR  RANDOM  NUMBER  GENERATION 
DIMENSION  ISIARI  (5)  ,FKS  (400)  ,  IFKS  (4 00 ) 

DATA  I  START  /7 184 1 57 45 , 1 039450965, 9 9 8356 979 , 
*124005191, 1035698051/ 

DATA  MULT  /32771/ 

DATA  ISIZE  /1 073741 824/ 

DATA  DNORM  /214748364 8. DO/ 

1  READ  (5,200)  Ml ,N 

READ  (5, 205)  VEL ,A ,TA, B,TB,C ,TC, M2 
READ  (5,215)  (STEP  (I)  , FKS  (I)  ,1=1 , N) 

C 

READ  (4 , 230)  KSET 
IRAND=ISIART (KSEI  +  1) 

C 

DO  50  IKS=1,N 

C  COMPUTE  RANDOM  NUMBER  BETWEEN  -.5  AND  +.5 
IRAND=MOD (IRAND*MULT , ISIZE) 

DR AND=DFLO AT (IRAND) /DNORM 
FN=FKS  (IKS) 

C  IF  MEAN  IS  ZERO  NO  VALUE  OTHER  THAN  ZERO  IS  POSSIBLE 
IF  (FN.LE.0.0D0)  GO  TO  50 
C  IF  THE  MEAN,  FN,  IS  GREATER  THAN  GAUSSL 
C  IT  IS  CONSIDERED  GAUSSIAN 

IF  (FN.GT. GAUSSL)  GO  TO  40 

C  USE  TRUE  POISSON  DISTRIBUTION  BY  CUMULATIVE  SUMS 
DN=DBLE  (FN) 

IF  (DRAND.LT. 0.0D0)  DR AND=DRAND+1 .0D0 
DP=DEXP (-DN) 

DSUM=DP 

C  THE  COUNT  USED  IS  THE  INDEX  WHEN  THE  SUM  EXCEEDS  RAND 
DO  20  1=1,100 

IF  (DSUM.LT. DfiAND)  GO  TO  15 
FKS  (IKS)  =FLO AT  (1-1) 

GO  TO  50 

15  DP=DP*DN/DFLGAI (I) 

DSUM=DSUM+DP 
20  CONTINUE 

WRITE  (6,100) 

100  FORMAT (•  FUZZ  -  UPPER  LIMIT  EXCEEDED') 

GO  TO  50 

C  MEAN  IS  TOO  LARGE  TO  USE  POISSON.  GAUSSIAN  IS  ASSUMED 
C  FORMULA  OBTAINED  FROM : 

C  APPROXIMATIONS  FOR  DIGITAL  COMPUTERS,  CECIL  HASTINGS  JR. 

C  PRINCETON  UNIVERSITY  PRESS,  PRINCETON,  NEW  JERSEY,  U.S.A. 
40  RAND=SNGL  (DfiAND) 

ETA  =  SQRT  (-ALOG  (RAND**2)  ) 


■ 


4 


iJ 


'£ 
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. 


X  =  (.01 0328*ETA+. 802853) *ETA+2.515517 

X=X/( (  ( .  00 1 308*ETA+ .  1892  69)  *ETA+ 1 . 432788)  *ETA+1 . 0) 

X=-SIGN  (EIA-X , BAND) 

FKS  (IKS) =FN+X*SyfiI (FN) 

IFKS  (IKS)  =IFIX  (FKS  (IKS)  ) 

50  CONTINUE 

WRITE  (6,200)  M1,N 

WRITE  (6,210)  VEL, A ,TA ,B , TB,C, TC , M2 
WRITE  (6,220)  (STEP  (I)  ,IFKS  (I)  ,1  =  1 , N) 

KSET=KSET+ 1 
IF  (KSET.EQ. 5) KSEI  =  0 
WRITE  (4 , 230)  KSET 
IF(N.GT.I)  GOTO  1 
STOP 

200  FORMAT  (16X, 12, 52X, 12) 

205  FORMAT (23X, F5. 3, 1 X ,F5 . 0 , 1 X ,F7 .3 , IX ,F5 . 0 , 1 X ,F7 . 3 , 
*1X,F5.0,1X,F7.3,1X,I1) 

210  FORMAT (23X , F5 . 3 , 1 X,F 6 . 0,F 8 . 3 , F6 . 0, F8 . 3 ,F6 . 0 , F8. 3 , II) 
215  FORMAT (  (5  (I3,F6.0,5X) ) ) 

220  FORMAT (  (5  (13 ,16,5 X) ) ) 

230  FORMAT  (12) 

END 


t 


. 


